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ABSTRACT 


This document is primarily concerned with the performance 
of the 13,3 GHz airborne scatterometer system which is used 
as one of several Johnson Space Center (JSC) microwave remote 
sensors to detect moisture content of soil. The system is an~ 
alyzed with respect to its antenna pattern? the signal flow in 
the receiver data channels, and the errors in the signal out- 
puts. The operational principle and the sensitivity of the sys- 
tem, as well as data handling are also described. No discussion 
is made about the interaction of the transmitted radio frequency 
(RF) waves with the terrain surface. The dielectric property 
of the terrain surface, as far as the scatterometer is concerned, 
is contained in the assumed forms of the functional dependence 
of the backseat taring coefficient of the incident angle 9 . 

The consequence of a minor difference in the methods of antenna 
pattern measurements between 1.6 GHz and 13.3 GHz scatterometer 
systems was pointed out in the document. It is also found that, 
because of the rapid antenna gain changes with incident angles 
in the 13.3 GHz system, the corrections to the variations of 
the aircraft parameters are crucial in data reduction. Fianlly, 
a few recommendations are made to improve the overall performance 
of the scatterometer system. 


iii 



ACKNOWLEDGEMENT 


Special thanks are due to S. C. Reid, who provided all of the 
laboratory calibration data and many helpful discussions. 


. . .iv 



CONTENTS 


Section 

1. INTRODUCTION 1-1 

2. SYSTEM DESCRIPTION 2-1 

2.1 SYSTEM SPECIFICATIONS 2-1 

2.2 THE ANTENNAE 2-3 

2.3 PRINCIPLE OF OPERATION 2-4 

2.4 THE TAPE RECORDING SYSTEM 2-10 

2.5 RECEIVER SENSITIVITY 2-12 

2.6 THE SCATTEROMETER DATA HANDLING 2-16 

3. SYSTEM ANALYSIS 3-1 

3.1 DATA PROCESSING EQUATIONS 3-1 

3.2 RECEIVER SIGNAL ANALYSIS 3-7 

3.3 THE ROLLOFF FUNCTION 3-13 

3.4 THE MEASURED SYSTEM PARAMETERS 3-18 

3.5 THE FUNCTIONAL DEPENDENCE OF L-BAND on 6 . . . . 3-19 

4. THE ANTENNA PATTERN 4-1 

4.1 THE MAIN BEAM GAIN 4-1 

4.2 THE SIDE LOBE LEVEL 4-8 

4.3 THE ANTENNA CROSS POLARIZATION 4-11 

4.4 THE VARIATION OF ANTENNA GAIN FACTOR WITH 6 ." . . . 4-18 

5. NUMERICAL RESULTS AND DISCUSSION 5-1 

5.1 THE COMPUTER PROGRAM 5-1 

5.2 THE DERIVATION 5-2 

5.3 A SAMPLE POWER SPECTRUM 5-6 


V 



Section Page 

5.4.1 ANTENNA PATTERN MEASUREMENTS 5-11 

5.4.2 MEASUREMENTS OF CALIBRATE CONSTANT AND 

ROLLOFF FUNCTION 5-11 

f 

5.4.3 AIRCRAFT PARAMETERS 5-12 

5.4.4 ADDITIONAL ERROR SOURCES 5-16 

5.4.5 A BRIEF COMMENT 5-16 

6. CONCLUSIONS AND RECOMMENDATIONS 6-1 

7. REFERENCES .' 7-1 

Appendix 

A. THE CALIBRATION DATA A-1 

A.l THE CALIBRATE CONSTANT A-1 

A. 2 VARIABLE GAIN RESPONSE TEST A-2 

A. 3 THE ROLLOFF FUNCTION A-4 

A. 4 THE RECIEVER NOISE FIGURE A-4 

B. THE FERRITE MODULATOR B-1 

C. THE MIXER MATHEMATICAL MODEL C-1 

D. THE COMPUTER PROGRAM FLOW CHART AND LISTING D-1 


vi 



TABLES 


Table Page 

2-1 THE 13.3 GHz SCATTEROMETER SYSTEM CHARACTERISTICS. . 2-2 

2-2 SPECIFICATIONS OP THE TAPE RECORDER/REPRODUCER 

SYSTEM USED FOR JSME (WIDEBAND GROUP I) 2-11 

2- 3 THE 13.3 GHz SCATTEROMETER SIGNAL TO NOISE RATIO . . 2-15 

3- 1 A NUMERICAL EXAMPLE OF THE BACKSCATTERED SIGNAL 

LEVELS AT THE RECEIVER OUTPUT 3-17 

4- 1 THE RELATIVE POWER CONTRIBUTIONS TO THE TOTAL 

BACKSCATTERED SIGNAL (dB) 4-19 

5- 1 A COMPARISON OF 0 ° VALUES (dB) DERIVED FROM 

TWO DIFFERENT APPROACHES - FORE BEAM 5-4 

5-2 A COMPARISON OF a° VALUES (dB) DERIVED FROM 

TWO DIFFERENT APPROACHES - APT BEAM 5-5 

5-3 ERRORS INTRODUCED IN THE 0° ESTIMATES DUE 

TO AN AIRCRAFT ROLL OF 5° 5-13 

5-4 ERRORS INTRODUCED IN THE 0° ESTIMATES DUE 

TO AN AIRCRAFT PITCH OF 5° 5-14 

5-5 ERRORS INTRODUCED IN THE 0° ESTIMATES DUE 

TO AN AIRCRAFT VERTICAL SPEED OF 4 M/SEC 5-15 

A-1 THE RECEIVER OUTPUT RESPONSES TO THE VARIABLE 

GAIN AMPLIFIER SETTINGS A-3 

A- 2 THE MEASURED RESPONSES (ROLLOFF FUNCTIONS) IN dB 

FOR THE 13.3 GHz SCATTEROMETER RECEIVER OUTPUTS. . . A-5 

D-1 THE COMPUTER PROGRAM LISTING , . . D-7 


vii 



FIGURES 

Figure Page 

2-1 The coordinate systems used in the 1.6-GHz 

scatterometer antenna pattern measurements 2-5 

* 

2-2 13.3 GHz scatterometer block diagram 2-6 

2-3 The block diagram of the 13,3 GHz scatterometer 

real-time data processing 2-9 

2-4 Channel assignment of signal outputs for three 

scatterometers 2-13 

2- 5 Scatterometer data handling flow chart 2-17 

3- 1 The resolution cell as defined by Doppler frequency 

bandwidth and antenna beamwidth 3-3 

3-2 Geometry of a scatterometer ground resolution 

cell 3-6 

3-3 The 13.3 GHz scatterometer rolloff functions. . . . 3-15 

3- 4 The functional dependence of 2^ on 9 3-22 

4- 1 The one-way gain contour of the 13.3 GHz 

scatterometer vertical receive antenna 4-2 

4-2 The one-way gain contour of the 13.3 GHz 

scatterometer vertical transmit antenna 4-3 

4-3 The antenna gain variation with incident angle. . . 4-4 

4-4 The antenna beamwidth variations with 

incident angle 4-7 

4-5 The cross-track two-way antenna gain profile. . 4-9 

4-6 The cross-track two-way antenna gain profile. . . . 4-10 

4-7 The cross-polarized gain of the vertical transmit 

antenna 4-12 

4-8 The cross-polarized gain of the vertical transmit 

antenna 4-13 

4-9 The cross-polarized gain of the vertical receive 

antenna 4-15 

4-10 The cross-polarized gain of the vertical receive 

antenna 4-16 


viii 



Figure Page 

4- 11 The variation of the antenna gain factor with* • • • 4-20 

5- 1 A power spectrum plot of a 13.3 GHz scatterometer 

data set - fore beam, over water .......... 5-7 

5-2 A power spectrum plot of a 13.3 GHz scatterometer 

data set - aft beam, over water 5-8 

5-3 A vs. 6 plot derived from figures 5-1 and 5-2 . . 5-10 

B-1 Ferrite modulator single sideband power B-3 

C-1 Characteristics of a nonlinear resistive element . . C-3 

D-1 The geometry of an area element used in the 

numerical integration D-3 

D-2 The MAIN program flow chart D-5 

D-3 Subroutine GAAT flow chart D-6 


ix 



LIST OF ACRONYMS AND ABBREVIATIONS 


Acronym 

AC 

C 

CAD 

cm 

CW 

dB 

dBm 

DC 

DDC 

EOD 

ESD 

FM 

g 

GDSD 

GHz 

GSFC 

Hz 

IR 

IRIG-A 

JSC 

JSME 

KHZ 

m 


Alternating Current 
Centigrade 

Computation and Analysis Division 

Centimeter 

Continuous Wave 

deciBel 

deciBels referred to 1 milliwatt 

Direct Current 

Data Distribution Center 

Earth Observations Division 

Engineering Systems Division 

Frequency Modulation 

Gram 

Ground Data System Division 
GigaHertz 

Goddard Space Flight Center 

Hertz 

Infrared 

Inter-Range Instrumentation Group Format-A 

Johnson Space Center 

Joint Soil Moisture Experiment 

Degree Kelvin 

KiloHertz 

meter 


X 



Acronym 


MIz 

MegaHertz 

NASA 

National Aeronautics and Space. Administration 

PCM 

Pulse Coded Modulation 

Pf 

picofarad 

REC 

Record 

RF 

Radio Frequency 

RMS 

Root Mean Square 

SSBM 

Single Side-Band Modulator 

TAMU 

Te^tas AceK University 

UK 

University of Kansas 

VSWR 

Voltage Standing Wave Ratio 

W 

Vertical-Vertical 

XMIT 

Transmit 


xi 



1. INTRODUCTION 


The NASA/JSC (Johnson Space Center) 13.3 GHz single-polarized air- 
borne scatterometer system has been extensively studied and used 
for various applications in the past years (refs. 1, 2 , 3^ 4, 5, 
and 6) . Krishen (ref. 1) and Bradley (ref. 2) used the data ob- 
tained from this system over the ocean to estimate the backscatter- 
ing coefficient at various incident angles and discussed the 
relation between and the ocean wind speed. Rouse (ref. 3) and 
Parashar et. al. (ref. 4) have demonstrated the capability of the 
scatterometer system in the sea ice classification work. Dickey 
et. al. (ref. 5) showed an increase of ..5-7 dB in the backscattered 
signals as the scatterometer flew from the dry to irrigated parts 
of the same agricultural field. The system performance studies of 
the scatterometer based on data measured prior to 1970 were made 
by Bradley (ref. 2) and by Krishen et. al. (ref. 7) . A more recent 
analysis of the same system was done by Rosenkranz (ref. 8) . 

Recently, the 13.3 GHz scatterometer and another two scatterometer 
systems at the frequencies of 0.4 GHz and 1.6 GHz were installed 
in a NASA C130 aircraft as the additional microwave sensors for 

the Joint Soil Moisture Experiment (JSME) . The antenna ’gain 

> 

pattern for the 13.3 GHz system was re-measured at JSC in early 
1976. To make a quick assessment of the data quality for each 
mission, a real-time data processor has recently been added to the 
system also. With these new feature and measurements in mind, it 
is the purpose of this report to re-evaluate the performance of 
the system. In particular, it will be pointed out that there is 
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a slight difference in the methods of antenna pattern measurements 
made at JSC for the 13.3 GHz system and at the Physical Sciences 
Laboratory (PSL) of the New Mexico State University for the 0.4 GHz 
and 1 . 6 GHz systems . As a consequence of this small difference in 
the measurement approach, the data processing equation for the 13.3 
GHz system would be different from that for the other two systems, 
unless the measured antenna patterns for all three systems were 
brought into conformity by certain coordinate transformations. 

This document is organized into ci:: csctions. Section 2 deals with 
the system characteristics, the antenna and the system functional 

I 

descriptions, the procedures of data handling, and the expected 
receiver sensitivity at aircraft altitudes. Section 3 gives the 
scatterometer system analysis. In this section, the data processing 
equations to estimate are presented and a typical backscattered 
signal is traced through the receiver channels. Some coefficients 
appearing in the data processing equations are identified with the 
newly measured system parameters in the process. The antenna 
main beam gain contour, beamwidth variations with incident angles, 
as well as the side-lobe and cross-polarization gain levels are 
examined and studied in Section 4. Section 5 gives the mimerical 
results. These include the estimates of a® by two different data 
processing equations, the presentation of a sample backscattered 
signal power spectrum data from a recent functional check flight "■ 
over Lake Livingston, Texas, and an estimate on the precision of 
the scatterometer system. Finally, the major conclusions of the 
report and a few suggestions on the future work are contained in 
Section 6. 
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Four appendices are included at the end of the document to provide 
a detail list of the results from the laboratory measurements on 
the 13 » 3 GHz scatterometer system parameters/ the mixer mathemati- 
cal model/ the function of the ferrite modulator/ and a listing of 
the computer program for the estimation of a®. 
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2 . SYSTEM DESCRIPTION 


2.1 SYSTEM SPECIFICATIONS 

The NASA/JSC 13.3 GHz airborne scatteroweter is a continuous-wave 
Doppler radar system (ref. 9 ) , designed to measure the backscatter 
ing coefficient a° as a function of the angle of incidence 6 for 
various types of the earth's terrain. The scatterometer antenna 

I 

pattern is fan-shaped which covers the 6 range, along the aircraft 
flight path, of approximately ±5° to ±60° from the nadir direction 
The scatterometer receives the backscattered signal at all angles 
of incidence simultaneously. As a result of the aircraft's for- 
ward motion, different Doppler frequency shifts are introduced in 
the return signals from different 9. The power of the return 
signals from a ground resolution cell associated with a given 6 is 
finally retrieved by bandpass filtering at the corresponding 
Doppler frequency during data processing. 

The overall system specification and capabilities are given in 

Table 2-1. The specification of the real-time data processor is 

not included, because the unit is still under modification and 

testing. The scatterometer is operated only in the vertical- 

^ ' 

transmit and vertical-receive mode. 

In the next few subsections, the antennae, the operational prin- 
ciple, and the receiver sensitivity are further elaborated. The 
tape recording system as well as the data handling technique are 
also briefly described. 
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TABLE 2“1 

THE 13.3 GHz SCATTEROMETER SYSTEM CHARACTERISTICS 


Type of System 

Transmit Frequency 

Transmitter Power Output 

Number of Antennae 

Antenna Gain (two-way) 

Antenna Beamwidth (two-way) 

Antenna Along-Track Coverage 

Transmitter/Receiver Isolation 

Antenna VSWR (Max . ) 

One-Way Side Lobe Level 

Maximum Skew Angle 
(Antenna Pattern) 

Polarization 

Dynamic Range 

Receiver Sensitivity 

Receiver Noise Figure 

Signal to Noise Figure 
(1500 ft. altitude) 


CW Doppler 
13.3 GHz ± 100 MHz 
1.0 Watt 

Two, one transmit and one receive 

~38 dB 

~3® 

-,120° fan beam 
~52 dB 
1.5 = 1 

<20 dB below main beam 
~±0.5° 

W 

~G5 dB 
-150 dBm/Hz 
20 dB @ 10 KHz 
>10 dB 
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2 . 2 THE ANTENNAE 


There are two identical antennae in the 13.3 GHz scatterometer 
system, one of them associated with the transmitter and the other 
with the receiver. Each antenna is made of a piece of slotted 
rectangular waveguide. Both antennae were installed underneath the 
aircraft body between the two wings, with the axes of the wave- 
guides perpendicular to the aircraft roll axis so that the antenna 
beams were shaped in the cross-track direction. The nominal one- 
way, gross-track beamwidth is about 3 to 4 degrees as measured from 
each of the uwo antenna radiation patterns. The along-track 
antenna beam coverage is normally from -60° in the aft to +60° in : 
the fore direction. There is no apparent aircraft structure 
present in the region from nadir out to the horizon in both fore 
and aft directions. Thus,' the data collected within both fore and 
aft antenna beam coverages should be of useable quality. 

Some difference in the methods of antenna gain measurements was 
noted for the 13.3 GHz and 1.6 GHz scatterometers . The coordinate 
systems used in these measurements were shown in Figure 2-1. The 
antenna gain measurements for the 1 . 6 GHz system were made with a 

conical cut' around the axis parallel to the flight direction (i.e., 

o * 

the angle 90 - 6 remains constant for a given cut) . / The same 

measurements for the 13.3 GHz, on the other hand, were made with 

a conical' cut around the axis perpendicular to the flight direction 

i.e., the angle 90° - remains constant for a given cut). Because 

of this difference, the data processing equations for the two 
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scatterometer systems would be slightly different/ as pointed out 
in Section 3. The antenna gain patterns for the 13.3 GHz system 
were measured with a mockup to simulate the aircraft frame. 

Throughout this report the incident angle 0 is always measured 
from the nadir, positive in the fore direction and negative in the 
aft direction. The cross track angle <|) is measured at the location 
of the antennae from the axis of the aircraft flight as shown in 
Figure 2-1 for the 13.3 GHz system. 

2.3 PRINCIPLE OF OPERATION 

The operational principle of the 13.3 GHz scatterometer had 
previously been described in several reports (refs. 2, 7, and 9). 

For the sake of signal analysis in the later sections and the 
completeness of this document, it will be briefly described again 
in this subsection. Figure 2-2 shows the block diagram of the 
scatterometer system. It is copied from the report of-Reid (ref. 9) . 

The scatterometer system continuously transmits and receives RF 
energy at 13.3 GHz through the transmit and receive antennae which 
respond strongly to the vertically-polarized waves only. The 
responses to the cross-polarized waves for both antennae are finite 
but negligibly small (see Section 4). Therefore, unless otherwise 
specified, the discussion throughout this document refers to the 
W (vertically-polarized wave transmission and reception) polari- 
zation state only. 
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Figure 2-1. — The coordinate systems used in the 1.6-GHz and 
13.3 GHz scatterometer antenna pattern measurements. 
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Figure 2-2.— 13.3 GHz scatterometer block diagram 




















The backscattered signals are received at all angles of incidence 
(fore and aft) simultaneously. The received signal is then split 
into two components by a hybrid coupler, one without phase shift 
(channel 1) and one with a 90° phase shift (channel 2) . A trans- 
mitter sample signal (local oscillator signal) is also fed to each 
of the two data channels through a directional coupler. The com- 
bined signals in channel 1 are amplitude modulated after passing 
through a ferrite modulator driven by an audio oscillator. The 
modulated signals are mixed in the mixer and the resultant output 
from this mixing is a zero IF frequency spectrum data whi ch i s 
amplified by a series of audioamplifiers and then recorded on a 
tape recorder. Except the amplitude modulation by a ferrite modu- 
lator, the signal in channel 2 follows the same steps. The resul- 
tant output signal is again a zero IP frequency spectrum data 
which is amplified by a series of audioamplifiers and recorder by 
the same recorder. The signals from the two data channels are 
recorded with the same head on the adjacent channels of the magne- 
tic tape to minimize the timing inaccuracy and skew effect. In 
either one of the two data channels, the fore and aft signals are 
folded in the mixer. As a result, the output signals from the two 
data channels are in phase quadrature with each other. During data 
reduction, these phase quadrature data are digitized and Fast 
Fourier transformed. Then the suiranation and difference of the trans- 
formed signal spectra from the two data channels are taken to re- 
cover the fore and aft Doppler signals. ORIGINAL I* AG E IS 

OP- POOE QUAUT? 

As indicated in Figure 2-2, the power level of the transmitter sample 

signal in Channels 1 and 2 were measured to be 37.4 dB and 36.5 dB, 

respectively, below the power level of the transmitter. The back- 

scattered signal in the receiver was also not split evenly 
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into the two data channels# being 3.6 dB and 3.3 dB in Channels 1 
and 2 respectively below the backscattered signal. The slight 
unevenness of both backscattered and transmitter sample signals 
would not affect the quality of the data, because the final gain 
adjustment could be made just before the data was recorded on the 
tape recorder. 

With the addition of the real-time data processor, the data quality 
of a given flight mission could be assessed quickly. Including 
this addition, the block diagram of the entire 13.3 GHz scatterometer 
system is shown in Figure 2-3. A junction box is provided to handle 
all the data flow. The backscattered signals in Channels 1 and 2 ' 
and the ADAS data, after passing through the junction box, are 
recorded on the tape recorder No. 1. This tape recorder also 
records data from both 0.4 GHz and 1.6 GHz scatterometer s as 
discussed in the next subsection. The backscattered signals and 
ADAS data are also channeled to the real-time data processor for 
immediate processing. The results from this processing are sent 
back to the junction box and plotted on a strip chart recorder. The 
4-channel range expander provides adequate dynamic range to handle 
the processed data. The ADAS and the processed data are also 
recorded by the tape recorder No. 2. 

To aid in the association of the backscattered signal at a given 
time to a ground location, a Zeiss camera system was installed in 
the aircraft to take the IR pictures of the terrain at approximately 
4 sec intervals. The time of the data and picture taking enters 
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Figure 2 - 3 .— The block diagram of the 13.3 GHz scatterometer real-time 


data processing. 



in both the tape recording system and the picture frames. By 
comparing the times recorded in both the magnetic tape and the IR 
pictures, the backscattered signal at a given Doppler frequency 
(or incident angle) can be identified with a ground location, 
subject to a reasonable uncertainty. 

2.4 THE TAPE RECORDING SYSTEM 

The output of the scatterometer is recorded by a Mincora Series 110 
Recorder/Reproducer System, Wideband Group I. The system is a 
solid state, compact, portable unic whi^^h hat. a bi-directional tape 
transport and directly interchangeable signal electronics plug-in 
modules for FM recording. It records data on a 1-inch magnetic 
tape with 14 data channels and with a carrier modulating frequency 
of 108 KHz. It has a bandwidth coverage from DC to 20 KHz. The 
signal to noise ratio is ~40 dB and the speed of the data recording 
is 30 inches per second. The specification of the system is given 
in Table 2-2. 

The input level of the tape recorder is set at 1 volt RMS nominal, 
but is adjustable to produce full scale modulation between 0.5 and 
10 volts peak to peak. The output level is 1 volt RMS nominal 
into a 50-ohm load. It is also adjustable from 0 to- 1.5 volts RMS. 
The input impedance is 10,000 ohms nominal shunted by 330 pico- 
farads maximum unbalanced to ground. 

The tape recorder No. 1 (see Figure 2-3) is shared by three air- 
borne scatterometers at the frequencies of 0.4 GHz, 1.6 GHz, and 
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TABLE 2-2. - SPECIFICATIONS OF THE TAPE RECORDER/REPRODUCER 
SYSTEM USED FOR JSME (WIDEBAND GROUP I) 


Magnetic Tape Size 

Tape Speed 

Tape Speed Accuracy 

Time Base Error 

Start and Stop Time 

% Flutter (Maximum 
Cumulative) at 30 ips 

Center Carrier Frequency 

Frequency Response 

Signal to Noise Ratio 

Input Level 

Output Level 

Input Impedance 

Output Impedance 
Ambient Temperature Range 
Thermal Drift 


DC Linearity 


Total Harmonic Distortion 


1" in width and 14" in diameter / 
containing 9,200' of tape 

30 ips (inches per second) 

0.03 ips 

0.6 sec at 50 KHz reference frequency 
5 sec. 

0.30 peak to peak at flutter band- 
width of 0.2 to 5 KHz 

108 KHz 

DC to 20 KHz 

40 dB (RMS/RMS) 

1 volt RMS nominal. Adjustable for 
input levels between .5 & 10 volts 
peak-to-peak 

1 volt RMS into a 50 ohm load. Ad- 
justable from 0 to 1.5 volts RMS 

10,000 ohms nominal shunted by 330 
pf maximum unbalanced to ground 

50 ohms maximum, unbalanced to ground 

0° to 50°C at sea level 

After 15 rain, warmup, less than 
±0.1%/C° in full scale output voltage 
change, over a temperature range of 
0° to 50°C 

Within ±0.5% of peak-to-peak devia- 
tion from the best zero based 
straight line 

< 1% at 1 KHz 


Humidity 


Up to 95% relative humidity without 
condensation 
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13.3 GHz. The assignments of the channel number for various data 
outputs of the scatterometers are shown in Figure 2-4. The 13.3 GHz 
scatterometer outputs are recorded on two adjacent data channels 
(11 and 13} of the same head. 

The time of data taking is recorded through IRIG Standard Time Code 
and Aircraft Data Acquisition System (ADAS) in channels 9 and 2 
respectively. In addition to the time information, ADAS also 
acquires the aircraft parameters such as altitude, speed, and roll, 
pitch and drift angles. 

The time and aircraft parameters provided by ADAS are recorded 
sequentially in PCM format and frequency modulated at 225 KHz. All 
of this information is derived during data reduction and, one way 
or another, is used in the data processing for the derivation of 
backscattering coefficient cr°. 


2.5 RECEIVER SENSITIVITY 

One way to estimate the sensitivity of the scatterometer system is 
by examining the signal to noise ratio at the receiver output. The 
signal at the receiver output can be evaluated by using .the data 
processing equation described in Section 4 and the measured system 
loss factor. For the convenience of a numerical estimate, Eq. (3-3) 
is repeated and simplified here: 


(4tt) 


f v! 

J R‘* 


5 P. •« 2 G. 0 G 

r _ tX-^ to 1 


dA = 


(4ir) ^ R'^ 


( 2 - 1 ) 
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13.3 GHz 
SCATT. 


1.6 GHz 
SCAH. 


0.4 GHz 
SCATT. 



Figure 2-4.— Channel assignment of signal outputs for three scatterometers , 


HINICOM 110 TAPE RECORDER/REPRODUCER SYSTEM 



where a bar over a given symbol indicates the average value of that 
parameter over a ground resolution cell A. The numerical values 
of all the parameters in the above equation were evaluated and 
entered in Table 2-3 for the estimate of the receiver output signal 
level. The aircraft altitude and speed were assumed to be ~460m 
and ~80 m/sec, respectively. The incident angle was taken to be 
60° where low signal return was expected. The along-track ground 
cell length used in the area calculation was ~25m which, at 60° 
incident angle and 80 m/sec aircraft speed, corresponded to a 
bandwidth of -48 Hs. The cross-track two-way beamwidth as obtained 
from Figure 4-4 was -3.9°. The power loss from the transmitter to 
the transmit antenna was measured to be -0.5 dB (S. C. Reid, personal 
communication) . The u° value at 60° incident angle was derived 

I 

from Figure 3-4; the smooth water valve was used for low signal 
return. All of the remaining parameters appeared in Eq. (2-1) are 
self-explanatory . 

The sum (is dB) of the radiated power and the propagation constant 

gives the received power at the antenna terminal. Since there is 

no receive antenna feed line loss, the receiver signal output is -120.9 

dBm, The amplification factor and the rolloff function were not 

* *■ 

entered in this calculation because both of them applied to the 
signal and noise equally and cancelled out when signal to noise 
ratio was taken. 

For the calculation of the receiver noise level, the ambient 
temperature was assumed to be -290° K. The Boltzman's constant 
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TABLE 2-3. - THE 13.3 GHz SCATTEROMETER SIGNAL 

TO NOISE RATIO 


Output Of the Transmitter 
Transmit Antenna Feed Line 
RADIATED POWER 

(Wavelength 0.02255 m) 
G^G^ (Gain at 0 = 60°) 

a° (at 0 = 60°) 

(1/4 TV)^ 

(1/R) ^ (Range) 

A (Area Element) 
PROPAGATION CONSTANT 
RECEIVER SIGNAL OUTPUT 
kT 

Bandwidth Corresponding to 
Cell A 

Receiver Noise Figure 
RECEIVER NOISE POWER 
SIGNAL TO NOISE RATIO 


+30 dBm 
Loss -0.5 dB 

+29.5 dBm 
-32.9 dB 
+27.0 dB 

-24.8 dB 
-33.0 dB 
-118.6 dB 
31.9 dB 
-150.4 dB 
-120.9 dBm 
-174.0 dBm 

Resolution 16.8 dB 

18.0 dB 
-139.2 dBm 
18.3 dB . ' 
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k is 1.38 X 10 joule/ K. The noise figure of 18 dB was obtained 
from the laboratory measurements at 6 KHz {see Appendix A) . Summing 
all of these factors, the receiver noise output over the 48 Hz 
bandwidth is estimated to be -139.2 dBm, which results in a signal 
to noise ratio of -18.3 dB. 

For operation over land or for small angle of incidence, the signal 
to noise ratio is expected to improve. For example, over land 
at 60® incident angle could be ~5.4 dB (see Figure 3-4) instead of 
-24.8 dB for- water. The signal to noise ratio in that case would 
be ~37.7 dB rather than -18.3 dB. The signal to noise ratio at 
smaller incidence angles could be better than 40 dB. However, 
the signal to noise ratio of the tape recorder system described 
in the previous subsection is -40 dB, which represents the best 
value attainable with the 13.3 GHz scatterometer and the data 
recording system. 


2.6 THE SCATTEROMETER DATA HANDLING 

A fair amount of data handling is required before the processed 
scatterometer data can be delivered to the principal investigators 
for studies and interpretations. Figure 2-5 gives a flow chart 
of the scatterometer data handling. The analog tape,' which con- 
tains the data for all three scatterometer systems (at frequencies 
of 0.4 GHz, 1.6 GHz, and 13.3 GHz) from a given mission, is 
duplicated at Ground Data System Division (GDSD) . The aircraft 
parameters such as roll, drift, pitch, altitude, and speed are 
also extracted and prepared in a tabular form at GDSD. The copies 
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td ^ Figure 2-5.““ Scatterometer data handling flow chart. 
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of the analog tape and the tables of aircraft parameters are 
delivered to the Engineering Systems Division (ESD) and the Data 
Distribution Center (DDC) , respectively for a quick-look processing 
and for distribution. The results from the quick-look analog 

* I 

processing are examined for data quality of the mission at the 
SF5 branch of the Earth Observations Division (E0D/SF5) . 

In addition to the recording of the scatterometer data and the 
aircraft parameters/ the flight logs and the Zeiss IR film are 
also generated in a mission. The flight logs contain information 
on the occurrence of events for the entire mission/ V7hile the IR 
film shows the ground scene of all the aircraft flight lines of 
interest. The IR film is processed and copies made at the JSC Photo 
Laboratory/ and the flight logs are duplicated at ESD. . Copies of 
both the IR film and the flight logs are delivered to DDC for 
further distribution. ESD also receives a copy of the IR film 
to aid in the quick-look data processing. 

Generally/ a system calibration is performed on the scatterometer s 

in ESD before a major mission takes place. The calibrate constants 

are needed in both quick-look and final data processings at ESD 

* ^ 

and Texas A & M University (TAMU) , respectively. TAMU receives a 
copy each of the analog tape, the aircraft parameters, the flight 
logs, and the calibrate constants to do the data processing with 
a digital method. The final outputs from the TAMU processor are 
plots and/or listings of the time history of the backscattering 
coefficient a°. These outputs are distributed to the principal 
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investigators at NASA/GSFC/ University of Kansas, and the University 
of Arkansas for further studies. E0D/SP5 also receives a copy for 
bookkeeping and evaluation purposes. 

The ground truth data collection is performed jointly by all the 
parties involved in the JSME mission. The collected data are 
processed and organized at the University of Kansas. The final 
ground truth data outputs are tabulated and distributed to all the 
principal investigators involved and to E0D/SF5 at JSC. 

j 

Finally, the flow and distribution of the 13.3 GHz scatterometer 
data processed by the real-time processor have not been determined 
by the time of this writing. The handling of those data may be 
decided in the near future. 
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3. SYSTEM ANALYSIS 


data processing equations 

power Pj;- reflected from a target and collected by a receiving 
can be expressed in terms of the radar parameters as 
. 12 ) : 

o 

^r (4TT) 3 R** (3-1) 

= power transmitted by the radar 

= gain of the transmitting antenna in the 
direction of the target 

Gj. = gain of the receiving antenna in the 
direction of the target 

a ~ radar cross-section of the target 

X = radar wavelength 

R = distance between the radar and the target 


? remote sensing of a surface terrain by an airborne radar 
-j^ystem such as the NASA 13.3 GHz scatterometer , it is more 
appropriate to use the backscattering coefficient (differential 
■'^xrr^ss section per unit area) rather than c. and o are related 

.^r.rough the differential area element dA by 

da = a° dA '(3-2) 

terms of this definition Eq. (3-1) becomes 




Pt Gt G,.j° 
(47T) 3 R'^ 


(3-3) 


3-1 



The power observed at a given central Doppler frequency 

and within a bandwidth Afa(0) is xaniquely related to a 

y .1 

ground cell A, as shown in Figure 3-1, through the relation 


f^(e) = ^ sin 0 

A 


(3-4) 


and 

Af^(0) = ^ COS0A0 (3-5) 

where V is the speed of the aircraft. The two isodopplers 
indicated in Figure 3-1 are defined by the two Doppler frequencies 
at f^j + Afjj and f^j - Af^. For numerical computations described in 
section 5, the ground cell A(0) is divided into N elements and 
Eq, (3-3) can be rewritten as 


Pj. (0) 



G^^(0) 


G^i(0) Ci' 


( 0 ) 


Ri^(0) 


aa^ 


(3-6) 


The power density Wj-(0) is simply 


Pr(Q) 

Wj.(0) = ^ 




N 


Af^(0) Af^(0) (4 tt) 


Z 

i = 1 


Gti(0) Gj-i(0) Ci°(0) 

R,‘*(0) 


Aa. (3-7) 


Two factors have been omxtted in the above discussion'. They 
should be considered before the final receiver output is obtained 
and interpreted. First, the attenuation of the RF power due to 
the two-way path of propagation is not included. The magnitude 
of this attenuation, in absence of heavy precipitation, depends 
on the amounts of oxygen and water vapor along the propagation 
path. At Skylab altitudes (''^235 nautical miles) and assuming a 


3-2 




Figure The resolution cell as defined by Doppler frequency 

bandwidth and antenna beamwidth. 





water vapor density of 7.5 gin/mS the two-way attenuation of the 
S-193 radar at 13.9 GHz and at incident angle was estimated 

to be only 0.24 dB (ref. 13). For most of the mission carried 
out by the C130 aircraft so far, the cruising altitude was 
•v 1,500 feet. Therefore, the loss due to atmospheric attenuation 
can be completely neglected. Secondly, Pj. as given in Eq. (3-6) 
is the power of the back-scattered signal at the receiver input. 

As the signal passes through the receiver data channels, it 
experiences losses as well as gains due to the presence of many 
subsystems. These gains and losses are automatically taken care 
of by the presence of a calibrate signal at 12 KHz as explained 
in the next few subsections. 

To derive the backscattering coefficient a° as a function of 

incident angle 0 from Eg. (3-7) requires an iteration process. 

First, a functional dependence of on 0 has to be assumed. Then 

the power density Wj. ( 0 ) of the backscattered signal is calculated 

from Eq. (3-7) for a series of 0 and A 0 . The calculated Wr(9) 

is compared with the power spectrum derived from observations. If 

the power spectra from computation and from observation do not match 

well over the incident angles of interest, another form of d°( 0 ) 

\ 

is assumed and the procedure repeated. The final functional 
dependence of a° on 0 is ascertained when best match between the 
computational and the observational power spectra is obtained. 

A simple, direct, and most used approach of determining as a 
function of 0 is to first associate a ground resolution cell 
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A(0) with a single average value of o°(e), as shown in Figure 
3-2, and then express o°(6) in terms of all measurable parameters. 
From this figure. 


and 


d^ = RA<{) 



RA6 

cos 


6 


A(6) 


di dj 


/ R«^'A4)A6 
COS 0 


(3-8) 


Substituting the expression for a 6 in Eq. (3-5) into Eq. (3-8) , 
the area o£ the ground cell becomes 

A(6) = (3- 

2V COS^0 


and the backscattered power from this area is 


0 3 

G. (0) G (0) a (0) X A(^Afd 

p (0) = _E — £ E. (3-10) 

^ (4tt) 3 2Vh^ 


Expressed in terms of c°(0), Eq. (3-10) becomes 


(4ti) 3 2Vh% (0) 

a°(0) = (3-11) 

P G^(0) G (0) x3a(1. 


where the aircraft altitude h = R cos 0 and Eq. (3-7) was used. 


Both approaches for the determination of a° versus 8 relationship 
are attempted in section 5 where the numerical calculation is 
performed. A comparison is made on the results from both approaches. 


3-5 
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The backscattered signal in the receiver data channels is further 
analyzed and the backscattering coefficient are discussed in 
the following subsections. The antenna gain pattern, side lobe 
level and coupling effect are discussed in detail in section 4. 


3.2 RECEIVER SIGNAL ANALYSIS 

To develop a comprehensive mathematical model for the signal 
response of the scatterometer system requires a complete know- 
ledge on the characteristics of each of the many subsystems 
shown in the block diagram of Figure 2-3. These characteristics 
include the gains, bandwidths, and noise figures of all ampli- 
fiers, the gains and bandwidths of all filters, to name just a 
few. To make new measurements on all of these parameters may 
be costly and time consuming, while to develop a mathematical 
model based on old values of these parameters may not be approp- 
riate, especially when some of the components were replaced over 
the past years. Consequently, a different approach is attempted 
here to derive a practical transfer function of the whole system, 
with the aid of some new measurements on a few key system para- 
meters . 

V 

N 

The approach follows closely with those attempted for the 13.3 
GHz airborne scatterometer system by Krishen et al., (ref. 7) 
and by Bradley (ref. 2 and for the 1.6 GHz system by Wang (ref. 11). 
The backscattered signal is followed through all major subsystems. 
The final outputs from the receivers contain factors such as the 
gains of the amplifiers, the coefficients of the ferrite modulators, 
and the gains of mixers and filters, which are combined and 
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identified with the key system parameters recently measured in the 
laboratory (see Appendix A) . 

Assuming that the continuous wave at frequency at the vertical 
transmitter port is represented by 


V. = V costo t 
too 


(3-12) 


The voltage at the input of the vertical receiver is 


V = Cj: cos ( (uj + tu,)t + X- I + C cos / (to - to-)t + x^ I (3-13) 
rr yod t J a lod a/ 


where 

= amplitude of the received signal in fore 
direction, in volts 

C *= amplitude of the received signal in aft 
^ direction, in volts 

to- = the Doppler shift corresponding to fore and 
aft beam at a particular angle of incidence, 
in radians/sec 

x^ = the phase shift introduced by the path length 
and the rough surface in the fore beam, in 
radians 

X = the phase shift introduced by the path length 
and the rough surface in the aft beam, in 
radians 


The direct leakage of the signal from transmitter to receiver at 
frequency to^ is not included in Eq. (3~i3) for reason given in 
section 4 . 
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The hybrid coupler in the receiver channel divides the input 
signal equally into two parts, one without phase shift and the 
other with a 90° phase shift. The output signal from the hybrid 
coupler can be written as 


cos ( (Wq + cos ( (Wq - ^a^ 


sin {(03^ + w^)t + x^) + sin ( (w^ - w^)t + x^) 


(3~14b) 


where and ^2c output signal voltages for the channel 

(channel 1) without phase shift and the channel (channel 2) with 
90^ phase shift, respectively, and k^ take care of the loss 

in the possible unbalance in the signal power division of the 
hybrid coupler. There exists an asymmetry between channels 1 and 
2 as shown in Figure 2-2. In channel 2, the backscattered signal 
from the hybrid output is directly mixed with the transmitter 
sample signal in a mixer. In channel 1, however, the transmitter 
sample signal and the backscattered signal are modulated by a 
12 KHz calibrate signal through a ferrite modulator. The signal 
from the ferrite modulator output is then mixed in the channel 1 
mixer. The voltage of the transmitter sample signal is 


= pV cos w^t (3-15) 

s o o 

and the signal voltage at the mixer inputs can be written as (see 
Appendix B) ; 
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ak- r 

Vjf Cj OOS ((<0^ + 0l^)t + Xf) + cos ((O^ - M^)t + Xg) 


. /5-PiV 

+ cos (D^t + ^ ; 


, cos o)_t cos o)^t I (3“16a) 

c o' 


^1 


ak. 


'2f 


r. 


= 7 =.|c^ Sin ((0)^ + (0^)t + Xf) + sin ( {(o^ - o)^)t + x^) 


>/^L 


0 '"d' 


/^P2^o 


COS (1) 




(3-16b) 


Here and P 2 gives the magnitudes of the amplitude attenuation 
in the transmitter signal for channels 1 and 2 respectively, 
is the angular frequency of the calibrate signal, a and ay (y<<D 
are the DC and AC gains of the ferrite modulator, respectively. 


After the mixers, the output signal voltages from both channels 
are (see Appendix C) : 


^lm = 




OJ^KtINAL PAGP Tg 

CfCos{m^t + x^) + C^cosCm^t - x^) op POOU 


2 YP^V 

+ ? cos (I) t 


V2m" 


1^2 P 2 


""o r . • 

— “CjSi 


sin(o),t + x<-) + C,sin(o),t - x ) 


] 


(3~17a) 


(3“17b) 


where 3^ and 32 take care of the conversion losses of the mixers 
in channels 1 and 2, respectively. The signal voltages as ex- 
pressed in Equations (3-17a) and (3-17b) are amplified by the 
series of audio-amplifiers in their respective receiver channels. 
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They also suffer attenuation through the presence of the switches 
and filters. If the total gains (before variable gain amplifiers) 
are and G^, respectively, the variable gain amplifiers are 
adjusted such that 


® “ ^l^lPl^o^l^lv 


^2^2^2^o®2*^2v 


(3-18) 


where G^^^ and 62 ^ are the gains of the variable gain amplifiers in 
channel 1 and 2, respectively. The signal voltages at the variable 
gain amplifier output are 



G 

/ 2 ~ 


^ YP^V 

C^COSCu^t + Xg) + C^cos(u^t - X^) + ^ cos UJ^t 

(3“19a) 



-Cj-sin(o 3 jt + X-;) + C sin(o),t - x ) 
r Cl jl a. cl a. 


(3-19b) 


The signals as expressed by Equations (3-19a) and (3~19b) are 
recorded by a tape recorder on two adjacent tracks of the same 
head. If additional gain unbalance is introduced by the tape 
recorder, it can be adjusted again in the data reduction. Since 
only the relative power levels of the backscattered and the cali- 
brate signals are used in the calculations of the backscattering 
coefficients, the gain adjustment does not affect the final results. 


The fore and aft Doppler signals are separated in the process of 
data reduction. The recorded data are first digitized by an 
analog-to~digital converter. The Fast Fourier Transform (PFT) 
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is then applied to the digitized data to convert the signals 
from the time domain to the frequency domain. The frequency 
spectra of the signals given by Equations (3-19a) and (3-19b) 
after the application of the FFT are; 


cos +, j6 (f - f^) 

V' [-CjU tf - ffl) 

sin + j6(f - f^) 


cos x^ - j6(f - f^) sin x^} + C^{6(f - f^) 


sxn X 




(3-20a) 


sin x^ + j5(f -■ f^) cos x^} + C^{-6 (f - f^) 


cos X,} 




(3-20b) 


Where j = \/-l, f ~ , and 6 (f - f^) is the Dirac delta function. 

'In the above expressions, the double-sided spectra (with both posi- 
tive and negative frequencies) after the FFT are converted to the 
single-sided spectra by the application of the power scaling factor 
of 2 (ref. 13). The signal voltage ‘densities , E, and E , (in volts 
per Hz) at Doppler frequency f^ in the fore and aft direction are 
obtained by forming the following expressions. 


s “ [(R^(Vi) + + (RjCV^) - 


(3-21a) 

(3-21b) 


The calibrate signal E (in volts per Hz) is given by 

c 


^c == 


Gyp^v^ 


(3-22) 
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From Equations (3-21b) and (3-22) , the received power density 
(in watts per Hz) is 


«ra<*d> - 4- 


= '(W (i) 


(3-23) 


Since the transmit power is proportional to , Eq. (3-23) can 
be rewritten as 


W, 


ra 




(3-24) 


where K is a measurable quantity in the laboratory, as discussed in 
subsection 3.4. 


3.3 THE ROLLOFF FUNCTION 

From Eq. (3-3) , it is seen that the power level P^. of the back- 

scattered signal depends on the backscattering coefficient a°, the 

distance R from the scatterometer to the surface area of interest, 

and the gains of the transmit and receive antennae. a° is generally 

observed to decrease with the angle of incidence 6 over either land 

or water. And R = h/cos6, h being the aircraft altitude, increases 

with 0. Since P„ is proportional to and inversely proportional 

to the fourth power of R, P^ is expected to decrease 'rapidly with 

0. Moreover, the antenna pattern studied in the next section shows 

a decrease in two-way gain with increasing 0 . This further enhance 

the decrease of P with O . The rapid decrease of P^ with 0 has an 

r r 

important bearing on the design and performance of the 13.3 GHz 
scatterometer system. For example, the' signal noise ratio of the 
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Minicom 110 Tape Recorder System described in subsection 2.5 is 
limited to a maximum of 40 dB. If at 0 = 5® is 40 dB or more 

higher than that at 0 = 60^, then the signal power level would be 
either saturated at e = 5° or buried in the noise at 9 = 60° during 
data recording. Therefore, it is necessary to install certain filters 
to offset the strong dependence of on 0. 

To reduce the power level of the signal return at small 0 before 
being recorded by the tape recorder, two rolloff filters were in- 
stalled in each of two receiver data channels. One of the filters 
was used for observations over lands and the other over waters. 

Either one of the two filters is made operational by a control 
switch. The responses of these filters as a function of frequency 
for both receiver channels were measured in the laboratory at 
Building 15 of JSC and given in Table A-2 of Appendix A. The 
frequencies at which the filter responses were measured can be 
converted to the angles of incidence through Doppler relation of 
Eq. (3-4) for a given aircraft speed. The listed numbers in Table 
A-2 are average values of measurements on both data channels at each 
of the frequencies indicated. These numbers were plotted as a 
function of frequency in Figure 3-3 for both land and water filters. 
The values at the frequencies corresponding to 0 » 5°, 10°, 15°, 

20°, 30°, 40°, 50° and 60° are obtained from these curves and used 
in the numerical calculation in Section 5. 

As an excimple in the application of the rolloff functions 2^(e> 
and ZA(0) given in Figure 3-3, observe that the values of a® over 
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calm water as obtained from Figure 3-4 in subsection 3.5 are _ 18 
dB and ^ -23 dB for 9 = 5® and 9 = 60°, respectively. Assming 
that the aircraft is cruising with a speed of _ 75 m/sec and at 
an altitude of ^ 480 m, the Doppler frequencies are 580 Hz and 
5760 Hz, and the distances are . 482 m and . 960 m for 9=5° and 
9 = 60°, respectively. The corresponding attenuations at these two 
frequencies are found from Figure 3-4 to be « 18 dB and 1.8 dB. 

If area A = 1 m^ , = 1 watt, and neglecting the signal losses 

and amplifications in the scatterometer system, the receiver out- 
put powers at 9 = 5° and 9 = 60° were calculated from Eq. (3-3) 

t, 

and given in Table 3-1. It is clear from the table that there is a 
48 dB difference in the signal returns at 0 = 5° and 9 = 60°, which 
is more than the dynamic range of a typical tape recorder. 

Therefore, on some mission over an extremely calm water surface, 
if the scatterometer variable gain amplifier is set to receive the 
backscattered signal return at 9 = 5° without saturation, then the 
signal return at 0 = 60° would be lower than the tape recorder noise. 

It was found out from many observations over lands that the use 

of water rolloff filters generally yielded better results for the 

13.3 GHz scatterometer. Therefore, it is reasonable ,to assume that 

the water rolloff filters will be most frequently used in the 

future missions either over land or over water. Only the water 

rolloff function Z (9) will be used in the numerical calculations 

(0 

in section 5. 
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TABLE 3-1. - A NUMERICAL EXAMPLE OF THE BACKSCATTERED 
SIGNAL LEVELS AT THE RECEIVER OUTPUT 


u> 

I 

H* 


^ w 

fir* 

8g 

Wfe 


Parameters 

0=5° 

0 = 60° 

Prom Eq. (3rl3) 

“ 70 Hz 

= 693 Hz 


1 Watt, ,30 dBM 

1 Watt, 30 dBM 


37 dB 

27 dB 

a° 

18 dB 

"-23 dB 


5x10“ V, -32.9 dB 

5xl0“V, -32.9 dB 

(4ir)3 

33 dB 

5.35 X 10^ V, 107.3 dB 

33 dB 

■8.46 X lO^V, 119.3 dB 

R** 

A 

Im^, 0 dB 

Im^, 0 dB 

P^{0) 

17 dB 
-105 dBM 

1.8 dB 
-153 dBM 



3.4 THE MEASURED SYSTEM PARAMETERS 


From the results of discussions in the previous three subsections, 
it is now possible to put all the parameters together and derive 
the final forms of the data processing equations. Note that both 
the’ backs cattered and the calibrate signals, E^ and E^ in Equations 
(3-21b) and (3-22) , pass through the water rolloff filters in the* 
receiver data channels and, therefore, subject to the attenuations 
of the filters amounting to Z^(f) (or Z^(e)). occurs at the 

frequency of 12 KHz so that the value of is ^ 0 dE, Conse- 
quently, only the Z^(f) for E^ needs to be taken into account and 
Eq. (3-24) should be rewritten as 


"ra«d> = 



(3-25) 


Combining Equations (3-6) , (3-7) , and (3-25) , the expression for 
the first form of the data processing equations becomes 




t' 1 


Gt^i) G^(0^) 


i = 1 




similarly, the second form of the data processing equation's is 
obtained by combining Equations (3-11) and (3-25) . The result is 


a'^(0) 


(4ti)^ 2Vh^Z (6) 
0) 

KGt(0)Gr(0)X^A4) 



(3-27) 
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Note that, in Eq. (3-*26) , Af^ is the relative receiver output 
power in the Doppler frequency band Af^ corresponding to the signal 
return from the ground cell defined by two isodopplers, f^ + 
and f^ - J£Af^. 

Two key system parameters were measured in Building 15 of JSC, 
namely, K, and Z^(6). The results of the recent measurements of 
these parameters were given in Appendix A. The rolloff filters 
were permanently installed in the scatterometer system. Conse- 
quently, the values of Z^(9), which were pre‘=!»=>nted in the previous 
section, are not expected to change much from measurements at 
different times unless the filters were replaced. The values of 
K, on the other hand, could change for different missions. The 
parameter K controls the power level of the calibrate signal which 
is also subject to the limitation of the tape recorder dynamic 
range. For a given surface of backscattering and aircraft alti- 
tude, K could be adjusted for an appropriate calibr ate signal 
level to bring about the best signal to noise ratio for the data. 
The values of K measured for Mission 347 functional check flight 
as presented in Appendix A will be used for the numerical calcula- 
tions in section 5. 

I 

3.5 THE FUNCTIONAL DEPENDENCE OF L-BAND g° ON 9 

It was pointed out previously that the ordinary approach of deriv- 
ing the backscattering coefficient a°, by means of a Doppler radar 
such as the 13.3 GHz. scatterometer, is to use the second form of 
the data processing equation given by Eq. (3-27) . The o® derived 
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from this approach contains two sources of error, both resulting 

from the approximation of a rectangular area (Figure 3“2) for an 

actual doppler cell defined by two isodopplers (Figure 3-1) - 

The first source of error comes from the assignment of a single 

value of for a given rectangular ground cell. Actually, 

was observed to vary with incident angle 9 , and the contribution 

to the observed signal return within a given bandwidth Af^ is 

the sum of the backscattered signals from all the small area 

elements with different e.'s in the doppler cell. 

^ I 

The second source of error is due to the crudeness of the area 
approximation itself. Both of the error sources depend on the 
scatterometer antenna beamwidth. It is desirable to explore the 
amount of error introduced in the derivation by using the 
second form of data processing equation. The derived from the 
first form of data processing equation given by Eq. (3-26) provide 
an adequate standard for comparison. 

o ‘ 

The derivation of a by Eq. (3-26) requires an assumed functional 

dependence of on 0 to begin with, as pointed out previously. 

To see how the shape of versus 0 curve affects the error in the 

\ 

estimation by Eq. (3-27) , the measured values of g® from 
both smooth and rough surfaces of backscattering are considered. 
For the smooth surface, the measurements made in February 1971 
over North Atlantic Ocean site 166 (mission 156) and reported by 
Krishen et. al. (ref. 7) are adopted. For the rough surface, the 
measurements made over a test site near Garden City, Kansas in 
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June 1970 (mission 133) and reported by Dickey et. al. (ref. 5) 
are used. The results from both of these measurements are shown in 
Figure 3~4. The land data shown in the figure (X sign) were de- 
rived from the smooth curve on Figure 4 of Dickey et. al. (ref. 5). 

The data over the calm sea showed a power law dependence of cf° on 
0. As displayed in Figure 3-4, the measured dependence on 6 was 
well represented by the relation; 

0 ° = 163.20“°*^® - 30.0 (3-28) 

A quadratic regression was performed on the data from the measure- 
raents over land. The result is given by: 

= 0.00270^ -0.35336 + 6.101 (3-29) 

In both of the above equations, the unit of a° is in dB and that 
of 0 in degree. Both of these expressions are used in the numeri- 
cal calculations in subsection 5.2, where the two forms of the 
data processing equations are compared. 
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Figure 3-4.— The functional dependence of 2** on 6* 
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4. THE ANTENNA PATTERN 


4.1 THE MAIN BEAM GAIN 

The most recent measurements of the gain patterns for both transmit 
and receive antennae of the 13.3 GHz scatterometer were carried 
out in 1976 in Bldg. 14 of JSC. The measurements were made for 
every 1° interval in both along-track and cross-track directions. 
Since it was known from the previous experience that the side-lobe 
gain (two-way) was ~30 dB or more below the main beam gain for these 
antennae, the gain patterns were measured only over ± 10° in the 
cross-track direction. The measurements were performed with mock- 
up to simulate the effect of the aircraft body. Consequently, the 
measured gain patterns are expected to be close to reality. 

To see the spatial variation of the antenna patterns, the contour 
maps of the one-way gains are displayed in Figure 4-1 and 4-2 for 
the transmit and receive antennae respectively. Each of the maps 
covers the angular ranges of <j) from ~-10° to ~+10° and of 0 from 
-70° to +70°. The direction of the nadir is located at 9=0 and 
(j) = 0. From these figures it is clear that both antennae show a 
narrow beamwidth of ~3° and three gain maxima within 0 range of 
+40°. The gain maxima of both antennae occur at 9 =.^-25°, -v0°, 
and ~+25°. Although the gain pattern measurements were limited to 
within ({) = ±10°, the data plotted in Figure 4-2 indicated a side- 
lobe level some 20 dB or more below that of the main beam. 

Figure 4-3 shows the variations of the antenna two-way gain with 
0 for 4* = Of “1° and +1°. The two-way gain maxima at 0 = 0° and 
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0~±25° are comparable at 36~39 dB. Additional gain maxima of 
smaller magnitudes (^31-35 dB) are also present at 6~±45°. Beyond 
0 = ±50°, the two-way gain decreases rapidly. Another feature dis- 
played by these curves is the asymmetry of the gain in the fore and 
the aft directions. For any 6 in the fore direction, the gain at 
(J) = 0° is smaller than those at ij) - ±1°. This implies a double 
main beam maxima in the cross-track gain profile which will be 
discussed in more detail in the next few subsections. 

To see the beamwidth variation with 9 , the same approach used in 
the analysis of the 1.6 GHz scatterometer system (ref. 11) is 
adopted here. Two methods of estimating the antenna beamwidth are 
defined. The first one is to take the maximum value of the two- 
way gain G^^F(6) at an incident angle 0. A 3-dB value was sub- 
tracted from G^^F(0) (in dB) . Two values of the cross-track angles 
and <^ 2 / measured from the location of G^^F(0), were determined 
where the antenna gain equaled G^^F(0) - 3 dB. The 3-dB beamwidth 
is the s\3m of <j>^ and <p 2 » The second ijfiethod (summation method) of 
defining a beamwidth is given by 

K"'G^(e,(t))G^(e,<}))d(j) 

A4)(0) = ^ (4-1) 

G^^F(0) 

where (j>^ and ^2 the cross-track angles where the null points 
of the main-beam antenna two-way gain occur. G, (9,()i) and G (6,^t) 
are the one-way gain at location defined by 0 and (}> of the transmit 
and receive antennae respectively. All the gains in Eq. (4-1) are 
expressed in numerical value (not in dB) . The results from those 
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two methods of bearawidth definition are shown in Figure 4-4. The 
dotted and solid curves in this figure give the beamwidth (A(J)) 
variation with 6 as determined by 3-dB and sxammation methods 
respectively. 

Clearly, there is an asymmetry in the variation of A(j) with 6 in the 
fore and aft directions. In the aft direction, A(j) decreases from 
.,4° at 0=~5° to ~2.4° at 0='-65*^. Throughout the 0 range from 
to -70°, the A^’s determined by these two methods track each other 
very well, the difference between them being not more than ~0.3°. 

On the other hand, A<{i's in the fore beam vary about a mean value of 
~4°. The range of variation is from 3° to 5° and from 2.1° to 5.7° 
for A^i determined by summation method and 3-dB method respectively. 
The reason for large variation in A(j> determined by 3-dB method is 
due to the double maxima in the cross-track gain profile mentioned 
previously. For large A^, the values of the two gain maxima differ 
by less than 3dB so that the locations of both maxima are contained 
within A<p . For small A<|) , the difference in the two maximum gains 
is more than 3 dB so that only the larger gain peak is contained 
in A(f). 

From the above discussion, it is clear that special care must be 
taken in the derivation of the backscattering coefficient from 
the data obtained in the fore beam. The bearawidth A<j> determined 
by 3-dB method could result in an appreciable error in the 
calculation. The A(J» determined by the summation method will be 
used in the numerical calculation presented in Section 5. 
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Figure 4-4.— The antenna beamwidth variations with incident angle 


4.2 THE SIDE LOBE LEVEL 


Since the antenna gain measurements were made only in the region 
within (J) - ±10^/ not many side lobes were observed in the measured 
data. Figures 4-5 and 4-6 showed, from the limited data, the 
cross-track gain profiles at two selected incident angles in aft 
and fore directions respectively. For 6 = -10° and 9 = -60° 
in Figure 4-5, the side lobe gains are down by ~55 dB and ~33 dB, 
respectively, compared to the main beam gains. In figure 4-6, the 
side lobe gain at 0 = 60° is again down by ~33 dB compared to the 
main beam gain. At 6 = 10°, there i» difficulty in determining 
the side lobe level. However, from the past experiences (S. C. Reid, 
personal communication) , the side-lobe gain at that incident angle 
is at least 50 dB below the main beam gain. The side lobe gain at 
other incident angles were also examined and were found to be at 
least -30 dB below the main beam gain maxima. 

It is noted that the antenna parameters G^(6), G^(0) and A^(6) 
appear together in the denominator of the data processing equation 
(3-27) . To give a crude estimate on the side-lobe contribution 
to the total input receiver signal, it is necessary to compare the 
relative values of the product G^^P^ (0) A(()^ (0) (antenna gain factor, 
subscript i refers to either side lobes or main beam) for the side 
lobes and their associated main beam. For example, the beamwidth 
of the main beam at 9 = -60° is -2.5° which together with the two- 
way gain plot shown in Figure 4-5 results in an antenna gain 
factor of -13 dB. The beamwidth of the first side lobe at 0 » -4° 
is -1.8° which combines with a two-way gain of --10 dB gives an 
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TWO-WAY ANTENNA GAIN, dB 



CROSS-TRACK ANGLE <^, DEGREES 

Figure 4—6. — The cross-track two-way antenna gain profile. 
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antenna gain factor of ~-25 dB. Assuming that the side lobes are 

I 

separated by and are of the same peak gain of ~~10 dB each, 

the total side-lobe contribution over the entire range of (j> from 
±5^ to ±90*^ is about -8.7 dB, which is -21.7 dB smaller than the 
main beam antenna gain factor . It can be shown that at other 0's, 
the total side-lobe contribution to the main beam signal is at 
least ~20 dB down and is totally negligible. 

Figure 4-5 also shows the more than 1^ difference in beamwidth 
between the small and large incident angles in the aft beam. And 
Figure 4-6 displays the double maxima in the cross-track gain pro- 
file in the fore beam as described in the previous subsection. 


4.3 THE ANTENNA CROSS POLARIZATION 

The antenna response to the cross-polarized RP waves at the same 
frequency of 13.3 GHz was not included in the most recent measure- 
ments of the gain pattern. The effect of the antenna cross- 
polarization response was found to introduce some error in the 

estimate of the cross-polarized backscattering coefficients , 

HV ! 

and for the 1.6 GHz scatterometer system; for the linear- 
polarized backscattering coefficients a® and c° , the effect was 
negligible. Since the 13.3 GHz scatterometer provides the estimate 
of cf^ only, it is expected that the error introduced by this 
effect should be small. To give a crude estimate of this effect, 
the measurements of the cross-polarized response made in 1969 were 
examined here. Figures 4-7 and 4-8 showed the measured cross- 
polarized gains as a function of 6 for the transmit antenna. 
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Figure 


■7.— The cross-polarized gain of the 
vertical tranmit antenna. 
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Figure 4-9 and 4-10 displayed the same type of data for the receive 
antenna. The measurements were made at a cross-track angular 
interval of 3°. The plots in Figures 4-7 and 4-8 showed that the 
high cross-polarization gains occured at 0=0*^, ±45° for all <j) and 
that higher gains were observed for (J>~0° than for 4>>Q°. The ob- 
served gain values were less than -10 dB. The cross-polarized gain 
plots for the receive antenna in Figures 4-9 and 4-10, on the 
other hand, showed the higher gain for 4>^0° than for (j><0°. All of 

the observed gain values were less than -15 dB. The prominent 

cross-polarized gain peaks occured at 0-0°, ±45° and ±55° for this 

antenna . 


The effect of the antenna cross-polarization response is to make 
undesirable contribution to the total receiver signal input. 
Including this cross-polarization effect, the receiver input power 
should be modified and expressed by 


^rvv ^ ® ^ 


X2p^AA(0) 

(4TT)3R'*(e) L 


G. (0)a° (0)G (0)+G. {0)a? (0)G (0) 

tvm' ' vv ' rvm' tvc ' hv rvm ' 


+ (6)u° (0)G (0)+G^ (0)0?, (0)G (0) 

tvm ' ' vh ' ^ rvc ' tvc ' ^ hh ' ' rvc ' ^ 


(4-2) 


where 


AA(0) = the area of the ground cell at incident angle 0. 


Gtvjj^O) = vertically-polarized gain of the transmit antenna 
at 0. 


4-14 



- - 6 ^ 



Figure 4-9,— The cross-polarized gain of the 
vertical receive antenna. 
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Figure 4-10.— The cross— polarized gain of the 
vertical receive antenna, 
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G^^^(8) = cross-polarized gain of the transmit antenna 


at 0. 


G^^(0) = vertically-polarized gain of the receive antenna 
at 6 . 


G^vc^®^ = cross-polarized gain of the receive antenna at 0. 

The remaining symbols in Eq. (4-2) were defined previously. The 

first term inside the bracket in the above equation is the dominant 

term which represents the backscattered signal of interest. The 

three remaining terms are contributions of the undesirable signals 

due to antenna cross-polarization effect. These three terms, in 

sequential order, represent: 

« The return signal at the receiver due to the signal 

transmitted in cross polarization by the transmit antenna 
and backscattered in vertical polarization. 


.• The return signal at the receiver due to the signal 
transmitted in vertical polarization by the transmit 
antenna and backscattered in horizontal polorization. 


® The return signal at the receiver due to the signal 

transmitted in cross polarization by the transmit antenna 
and backscattered in horizontal polarization. 

Based on Figures 4-7, '4-8, 4-9 and 4-10, and the main beam 
antenna gain patterns, the relative contributions to the total 
power return of four terms in Eq. (4-2) can be estimated at some 
selected locations. The results of these estimates at six chosen 
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locations, where prominent cross-polarized gains were observed, were 

given in Table 4.1. The estimates were made with respect to the 

product of the two-way gain and o°. It is clear from this table 

that, without considering the relative values of all the o°*s, the 

signal of interest as represented by G. o® G is about 20 dB or 

tvm vv rvm 

more greater than the three undesirable signal contributions due to 
antenna cross-polarization effect, is generally comparable to 

a^h least a few dB greater than and This tends 

to make the undesirable signal contributions even more insignificant. 
Although the estimates were only made at six locations, it is quite 
clear that, if the backscattered signals were integrated over a 
ground cell as was done in the 1.6 GHz scatterometer system study 
(ref. 11) , the contributions due to antenna cross-polarization 
response would be totally negligible. 

4.4 THS VARIATION OF ANTENNA GAIN FACTOR WITH 6 

It was noted previously that the antenna gain G^^F{0) and beam- 
width A(j>(6) appeared together in the denominator of Eq. (3-27). 

This combined antenna gain factor directly affects the estimate 
that its variation with 0 should be examined in more detail. Figure 
4-11 shows the dependence of this factor on 9 in both fore and aft 
directions. It is seen from this curve that there are six or 
seven ripples in the 0' range from -60° to +60°. The existence of 
these ripples would require accurate d<^ termination of 0 during 
data processing. For example, if the error in the pitch angle 
correction is ~1°, then the error in a° at some 9's could be as 
large as ~1 dB. 
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TABLE 4^1 


THE RELATIVE POWER CONTRIBUTIONS TO 


THE TOTAL BACKSCATTERED SIGNAL (dB; 


\ 

KO 


II 



Along-Track 
Angle 9 , 
Degrees 

Terms 

Cross-Track Angle 4)/ 

Degrees 

0 

+3 

-3 


G. G 0® 

tvm rvm vv 

C 29.7 

<v 29.7 

“w 29-8 

+45 

G 1 G 0, 

tvc rvm hv 

“hv - 

“hv - 


G. G 0°, 

tvm rvc vh 

<h - 


*• ^-2 


O C' ^ 

tvc rvcahh 

“to - 

“hh -”•! 

" 28.4 


G. G 0 

tvm rvm vv 

0° +34.7 

w 

21.9 

“w + “-1 

-45 

^tvc^rvm^hv 

“hv * 

<v “ "‘0 

0? - 5.5 

hv 


^ tvm*^r VC ^ vh 

“vh - 

<h - 

“vh - 21-7 

i 

s 

G. G 

tvc rvc hh 

“Sh - 

0°h - 42.0 

“hh - 



ANTENNA GAIN FACTOR 6^F(e)Ai^(e) ,db 



ANGLE OF INCIDENCE e. DEGREE 

Figure The variation of the antenna gain factor with 9. 



The rapid variation of G^^F(9} A(j) (6) with 9 also requires that the 
mounting of the antennae onto the aircraft should conform to the 
laboratory mounting when the gain patterns were measured. Another 
words, the coordinate system in the laboratory gain measurements 
should be precisely the same as the aircraft coordinate system. 

If imprecision exists in the antenna mounting either in the labora- 
tory or the aircraft, error would result in the o® estimate which, 
if large enough, would make the interpretation of the backscattering 
phenomena difficult. 
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5. NUMERICAL RESULTS AND DISCUSSION 


5.1 THE COMPUTER PROGRAM 

The numerical calculations of the receiver signal outputs were per- 
formed according to Eg. (3-26) . The area elements used in the 
numerical integration covered 1^ in The centers of the area 
elements lie on the isodoppler corresponding to the along-track 
incident angle 6 . The area element dAj and distance R. have to 
be expressed in terms of the measured 0 and <f) . The derivation of 
the expressions for these two parameters were given in Appendix D. 
The antenna gains 0^(6^^^) and G^(6^), the backscattering coefficient 
a®(0^) and the distance R{6^) within the summation sign of Eq. (3- 
26) were all evaluated at the centers of AA^'s. When these centers 
did not fall on the positions where the actual measurements of 
GfCSi) and ^^.(0^^) were made, linearly interpolated values from the 
immediate neighbors were used. The summation was carried out to 
the null points on both sides of the main beam. 

The computer program was designed and implemented in accordance 

with Eq. (3-26) . It evaluated the power outputs in terms of the 

square of the backscattered to calibrate signal voltage ratio, 

\2 

I . Eight angles of incidence at 5°, 10 , 15®, '20®, 30®, 40®, 
o 

, and 60 were chosen for the computations. Since the along- 

track cell length is a fixed value of 25m, the bandwidth Af^ 

changes with angle of incidence, decreasing from 381 Hz at 0 = 5® 

o ® 

to 48 Hz at 0 = 60 . Eight values of rolloff functions 
corresponding to the eighty incident angles are entered in the 
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program as a table. The value of the calibrate constant K was 
given in Appendix A and treated as a constant in the program. The 
effects of both the antenna cross-polarization response and the 
side lobe gain were already shown in Section 4 to be negligible 
in the receiver signal input. Thus, the computer program did not 
include the computations of these two undesirable signal contri- 
butions . 


The program was coded in Fortran. It consisted of the main program 

and a subroutine. All of the computations leading to the v*»?ues 
2 

of ( I for all selected incident angles were performed in the 


(%) 


the incident angle 0^ and cross-track angle 


main program. The subroutine was used to obtain the values of the 
two-way antenna gain for the calculation of signal return, given - 

The two-way 

antenna gain pattern was previsouly saved on a magnetic tape and 

/ 

read and stored in the disk during the program run. The program 
listing and flow chart are given in Appendix D. The program was 
run on the UNIVAC 1110 computer at JSC. 


5.2 THE a DERIVATION 


As pointed out in Section 3, it is necessary to assiame a^ functional 

f 

dependence of on the incident angle 0 in the first method of 

deviation. This assumption was made according to Egs. (3-28) and 
(3-29) for the backscattering from water and land, respectively. 


With this assxmption the receiver signal outputs in terms of 


) 


Ea '2 


could be obtained by running the program listed in Appendix D. The 


signal outputs obtained in this way could be compared with those 
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actually obtained from observations to determine the validity 
of the assumed dependence on 6 . The procedure is straightfor- 
ward and will not be performed here. Instead/ a comparison will 
be made between the first (as represented by Eq. (3-26) ) and second 
(as represented by Eq. (3-27)) methods of a® derivation. 

/®a 

To do such a comparison/ the values of ( g— I resulted from the 
computer runs were taken to be the receiver output signals for 
both methods. Then the a® dependence on 0 as derived from the 
first method will be given by Eq. (3-28) for water and Eq. (3-29) 
for land. The numerical values from these equations at the 
incident angles of ±5°, ±10°/ ±15°/ ±30°, ±40°, ±50°, and ±60° 
were listed under the "a°" columns in Tables 5-1 and 5-2. Table 
5-1 gave the 0 ° values for the fore beam and Table 5-2 for the aft 
beam. a°'s from both over water and over 'land were entered in 
these tables. 

The corresponding values of the a° ' s derived from the second 
method of Eq. (3-27) were given in the same two tables under the 
"O 2 " columns. The same numerical values for the parameters such 
as altitudes, aircraft speed, cell length, and calibrate constant 
K were used in Eqs. (3-26) and (3-27). It is clear from these 
tables that, for observations over land where a° vs 9 curve is not 
steep, the estimated from both methods are comparable for all 0 
selected for computation. For backscattering from smooth water 
surface, the estimated a°'s from the two methods could differ 
appreciably for small incident angles, as can be seen from these 
two tables. For |0l ^ 20° where 0 ° vs 0 curve (see Figure 3-4) is 
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TABLE 5-1. - A COMPARISON OF 0° VALUES (dB) DERIVED 
FROM TWO DIFFERENT APPROACHES - FORE BEAM 


Angle of 
Incidence 0 
(Degrees) 

Wat 

.er 

Land 



o 

^1 

“2° 

5 

18.03 

15.78 

4.40 

4.27 

10 

-1.64 

-2.08 

2.84 

2.78 

15 

-9.16 

-9.32 

1.41 

1.38 

20 

-13.25 

-13.33 

0.12 

0.09 

30 

-17.69 

-17.73 

-2.06 

-2.09 

40 

-20.11 

-20.12 

-3.71 

-3.72 

50 

-21.65 

-21.73 

-4.81 

-4.89 

60 

-22.73 

-22.76 

-5.38 

-5.40 











TABLE 5-2. - A COMPARISON OF cr° VALUES fdB DERIVED 
FROM TWO DIFFERENT APPROACHES - AFT BEAM 






not steep, the estimated values from the two methods are again 
comparable. While for |6| £ 10° where the a° vs 0 curve is steep, 
the second method based on Eq. (3-27) tends to give a lower c° 
value than the first method. At 0 = 5°, the difference in a° 
estimate^ from the two methods could be as much as ~2 dB. 

Prom these comparisons, it can be concluded that when the a° is 
steeply decreasing with 0, the second method of a° derivation tends 
to underestimate the a° values. This suggests that, for observa- 
tion over a smooth quiet water surface, the a° values at small 0 
derived from the second method (which is used in most cases) may 
have to be modified. 

5.3 A SAMPLE POWER SPECTROM 

One way to show the scatterometer receiver data output is to plot 
the signal power density as a function of frequency. Figures 5-1 
and 5-2 gave samples of such plots for the fore and aft beam 
respectively. The data shown in these figures were obtained from 
a functional check flight over Lake Livingston, Texas, in June 
1977. The average aircraft speed and altitude were 77 m/sec and 
469 m respectively. It was rather windy and the water surface was 
not smooth during the flight. The power spectra in both" of the 
figures were generated with a bandwidth of 50 Hz and an integration 
time of 6 sec. 

Several features are apparent from the observation of these figures. 
First, the power spectra from both fore and aft beams are comparable 
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in general shapes and magnitudes. The noise floor was about ~~47 
dB in both figures. Secondly, the fore-beam power spectrum showed 
several distinct ripples in the frequency range of interest <5 KHz. 
The aft~beam power spectrum in Figure 5-2, on the other hand, showed 
a reasonably smooth variation over the same frequency range. The 
difference in the fore - and the aft - beam power spectra could be 
due to the direction of the wind blow relative to the aircraft 
flight line. Thirdly, the signal clipping occurred at the frequency 
of ~150 Hz in both figures, although the water rolloff filters were 
used during the flight. Of course, the gain can be adjusted to 
avoid the signal clipping, but this can be done only at the expense 
of the small signal to noise ratio at large incident angles. The 
small dynamic range of the rolloff filters were previsouly noted 
in Section 3. 

The backscattering coefficient a® could be estimated from Figures 
5-1 and 5-2 at several selected incident angles 6 . This was done 
using Eq. (3-27) and the results were plotted in Figure 5-3. The 
data obtained by Krishen (ref. 1) using the same scatterometer 
system was also included in the figure for comparison. The a^’s 
derived from the fore and aft beams were comparable to within 
dB for most 0. They were also not too different from 'the data 
measured during a wind speed of ..25 knots (dashed curve) . The 
a® values at 0 = 50° might not be good because of the small signal 
to noise ratio. 
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5.4 PRECISION IM THE OUTPUT SIGNAL 


The discussion in this section is limited to those factors directly 
related to the 13.3 GHz scatterometer system only. The backscat- 
tered signal statistics which were discussed in great length else- 
where {ref. 12) will not be repeated here. The following gives a 
brief description for each of these factors considered to affect 
the performance of the scatterometer system. 

5.4.1 Antenna Pattern Measurements 

The uncertainty in the antenna pattern measurements was not recorded 
and supplied together with the gain values. A reasonable estimate 
placed this measurement uncertainty at -±0.5 dB (S. C. Reid, personal 
communication). An uncertainty of ~±0.5 dB in one'-way gain results 
in the uncertainty of ~±0.7 dB for the two-way gain used in the 
estimates. Another possible source of uncertainty comes from the 
mounting of the antennae onto the aircraft. The antenna gain pat- 
tern was measured in a JSC laboratory. If care was not taken in the 
mounting of the antennae in the laboratory or in the aircraft, the 
laboratory coordinate system might not coincide with the aircraft 
coordinate system. As a result, a constant mounting bias might 
exist in the angle of incidence at which the gain was use'd in the 
estimates. The uncertainty from this source can not be easily 
assessed. 

5.4.2 Measurements of Calibrate Constant and Rolloff Function 

An uncertainty also exists in the measurements of the calibrate 
constant K and rolloff function Z^(f^) . A reasonable estimate 
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placed this uncertainty to be no more than ±0,5 dB for each of 
these measurements (S. C. Reid, personal communication) , 

5.4.3 Aircraft Parameters 

Factors associated with aircraft motion such as roll, pitch, drift, 
vertical velocity, speed and altitude changes could introduce 
errors in the received signals. The effects of roll, pitch, and 
vertical velocity are to alter the incident angle 9. These effects 
were studied by Krishen et. al. (ref. 7) with respect to the same 

13.3 GHz scatterometer system using the antenna pattern data 
measured in 1969. Following the same approach, the errors intro- 
duced by these factors in the derivation were estimated with 
respect to the newly measured antenna pattern data. The results 
were presented in Tables 5-3, 5-4, and 5-5. in making these error 
calculations, the aircraft speed is assiimed to be 80 m/sec, the 
vertical velocity to be 4 m/sec, and both roll and pitch angles to 
be 5°. Clearly, the change in 0 due to pitch and vertical speed 
could be either plus or minus. The examples shown in Tables 5-4 
and 5-5 consider the case of larger modified incident angles in 
the fore beam only (smaller modified incident angles in the aft 
beam) . The roll angle also changes the distance R. The values 
shown in Table 5.3 also includes the effect of the distance change. 

The changes of both aircraft speed and drift would introduce 
changes in the Doppler frequency associated with a given 0. For 
example, at 0 = 30°, a change of speed from 80 ra/sec to 75 m/sec 
would give a corresponding change in Doppler frequency from 3548 Hz 
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TABLE 5-3. “ ERRORS INTRODUCED IN THE cf° ESTIMATES 
DUE TO AN AIRCRAFT ROLL OF 5° 


Fore Beam 

Aft Beam ■ 

Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Error in 
(dB) 

Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Error in 
a® 
(dB) 

5 

7.07 

•>0.9 

-5 

-7.07 

-1.9 

10 

11.17 

-0.2 

-10 

-11.17 

0.2 

15 

15,79 

0.6 

-15 

-15.79 

0.1 

20 

20.59 

0.2 

-20 

-20.59 

- 0.1 

30 

30.38 

-0.3 

-30 

-30.38 

-0.4 

40 

40.26 

0.1 

-40 

-40.26 

-0.2 

50 

50.18 

-0.2 

-50 

-50.18 

-0.3 

60 

60.13 

-0.1 

-60 

i 

-60.13 

-0.2 
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TABLE 5-4. - ERRORS INTRODUCED IN THE a° ESTIMATES 
DUE TO AN AIRCRAFT PITCH OF 5*^ 


j Fore Beam 

Aft Beam 

. Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Errog in 
(dB) 

Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Errog in 
(dB) 

5 

10 

-3.3 

“5 

0 

2.4 

10 

15 

0.8 

-10 

-5 

2.8 

15 

20 

2.0 

“15 

-10 

“2.4 

20 

25 

0.4 

-20 

-15 

-0.1 

30 

35 

“3.5 

“30 

-25 

3.1 

40 

45 

-0.1 

“40 

-35 

-1.5 

50 

55 

“4.0 

“50 

-45 

4.1 

60 

65 

“3.5 

“60 

-55 

2.4 







TABLE 5-5. - ERRORS INTRODUCED IN THE 0 ^ ESTIMATES 
DUE TO AN AIRCRAFT VERTICAL SPEED OF 4 M/SEC. 


Fore Beam 

Aft Beam 

Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Error in 
0 

(dB) 

Angle of 
Incidence 
(Degree) 

Modified 
Angle of 
Incidence 
(Degree) 

Error in 
0 

(dB) 

5 

7.86 

~ 1.1 

-5 

-7.86 

-2.0 

10 

12.86 

-0.3 

-10 

-12.86 

X. • X 

15 

17.86 

1.5 

-15 

-17.86 

0.1 

20 

22.86 

0.1 

-20 

-22.86 

0.8 

30 

32.86 

-1.9 

-30 

-32.86 

-2.1 

40 

42.86 

1.0 

-40 

-42.86 

-0.5 

50 

52.86 

-2.2 

-50 

-52.86 

-2.0 

60 

62.86 

-1.4 

-60 

-62.86 

-2.5 



to 3326 Hz. Similarly/ a drift of ~10° would shift the Doppler 
frequency from 3548 Hz to 3494 Hz. These changes in Doppler 
frequency, if not corrected for in data processing, could result in 
erroneous a® estimates. 

The variation in aircraft altitude would also introduce error in 
the estimate. From Eq. (3-3) , is observed to be proportional 

4 

to R and inversely proportional to the area A of the ground cell. 

2 o 

Since R is proportional to altitude h and A to h , varies with 
h^. Thus, if h changes from 460 m to 440 m, the error in the 
estimate would be -.0.39 dB. 

5.4.4 Additional Error Sources 

Several additional error sources should also be considered. The 
first source was discussed by Krishen et. al, (ref. 7) in great 
detail. The second source of error had its origin in the elec- 
tronics assembly. Zusfeldt (ref. 14) analyzed such an error source 
in detail for the 1.6 GHz scatterometer system. It is expected that 
the same kind of error also exists in the 13.3 GHz system. How- 
ever, the magnitude of this error source is small and is mostly 
taken care of by the system calibration in the laboratory. Finally, 

■V 

if vs. 6 curve is steep, error could be present in estimate 
at small incident angles, as pointed out in Subsection 5.1. 

5.4.5 A Brief Comment 

If the mounting of the antennae in the laboratory gain measurements 
or in the aircraft should not be precise, then error would be 
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present in the a° estimates. Because of the existence of many 
ripples in the antenna gain factor variation with incident angle 
y Bt this imprecise antenna mounting tends to create bumps in the 

/C 

c° vs. 6 curve. As a result, the interpretation of the data with 
respect to the wave backs cattering from terrain surface may turn 
out to be very difficult. Thus, the error in the antenna mounting 
should be kept as small as possible. 

The calibrate signal level could be adjusted to accomodate a par- 
ticular mission. For each adjustment a nevr calibrate constant K 
would have to be measured and the measurement, therefore, should 
be made as precise as possible. The changes in aircraft para- 
meters could introduce large error in estimation. This type 
of error could be minimized by making corrections to the aircraft 

parameters during data processing. The error sources from the 

o 

improper phase shift and the steepness of o vs. 0 curve could 
be reduced to a minimum by proper data reduction and manipulation. 
However, if the data quality is such tha't the ratio of the fore 
to aft signal powers is large or small compared to 1, a few degrees 
in phase shift error could result in large error in estimate. 

The data should be used with reservation in this case. 


ORIGINAL PAGE IS 
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6. CONCLUSIONS AND RECOMMENDATIONS 


A system analysis was performed for the JSC 13.3 GHz airborne 
scatterometer . The principle of the scatterometer operation, 
calibration procedures, and data handling were described in some 
detail. The antenna gain patterns were studied in great length 
and the effects of the antenna side lobes and cross polarization 
analyzed. Several error sources in the estimation of were 
briefly discussed. In the course of the entire analysis, only the 
immediate power input to and the power output from the scattero~ 
meter system were dealt with. The interacting process of the RP 
waves with terrain surface was completely omitted. The terrain - 
surface characteristics were reflected through the assumed func- 
tional forms of the backscattering coefficient a*^(0). 

The major conclusions that resulted from this study are listed in 
the following: 

1. When 0 ° decreases rapidly with incident angle 0, which 
occurs in the backscattering from a quiet water surface, 
the normal approach of estimating 0 ° appears to give a low 
value at a small incident angle. If 0 ® vs 0' curve over 
water is given by Figure 3-4, then the amount of the 
underestimate is ~2 dB at 0 = 5°. 

2. The contributions to the receiver signal input due to the 
effect of the antenna side lobe gain and cross polarization 
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were examined in detail. These undesirable signal 
contributions were found to be negligible. 

3. It was found that the rolloff filters of the 13.3 GHz 
scatteroraeter system did not cover enough dynamic range. 
Because of this lack in dynamic range, the signal power 
input at low frequencies (small incident angles) might 
saturate for observation over a smooth surface. If the 
gain was adjusted to avoid signal saturation, then the 
signal to noise ratio would deteriorate at high fre- 
quencies (large incident angles) . 

4. The method of antenna gain measurements for the 13.3 GHz 
scatterometer was found to be somewhat different from 
that for the 1 . 6 GHz scatterometer . As a result , the 
data processing equations for the two syst^s should 
differ also. To maintain the same processing equation 
for both systems, the measured antenna gain pattern of 
one system would have to be coordinate - transformed in 
conformation to the other. 

5. Because of the rapid changes of the antenna ga'in factor 
with the angle of incidence, the variations of the air- 
craft parameters were found to strongly affect the 
estimation of <j®. For example, a pitch of 5*^ could 
produce as much as ~4 dB difference in a® estimate at 
some incident angles. 
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6. The antenna gain pattern in the fore beam was found to 
possess double maxima in the cross-track gain profile. 
This would present some difficulty in determining the 
antenna beamwidth with the 3-dB method. 

There are five recommendations that are considered to be important 
in the improvements of the system performance and the quality of 
the final results. They are listed below: 

1. For observations over a smooth surface where a® decreases 
rapidly with 6 , care should be taken not to underestimate 
the a® values at small 6. 

2. The sea rolloff filters in the 13.3 GHz scatterometer 
system have been utilized in both land and sea observa- 
tions. It appears that the land rolloff filters are not 
of much use to the system. The sea rolloff filters alone 
do not have enough dynamic range to cover observations 
over both rough and smooth surfaces. Perhaps the land 
rolloff filters should be modified to improve the dynamic 
range of the entire system. 

3. The mounting of the antennae onto the aircraft or in the 
laboratory for gain measurements should be done with 
extreme care to minimize error. A mounting error of 1° 
could introduce ..1 dB error in a® estimates at some 
incident angles. Once the antennae were installed in the 


6-3 



aircraft, perhaps they should not be taken out again 
unless it was absolutely necessary to do so. 

4. Although the operational principle is the same for both 
the 13,3 GHz and the 1.6 GHz scat terome ter s, the slight 
differences in the laboratory calibration (see Appendix A) , 
in the antenna gain measurements, and in the location of 
the calibrate signals should be noted. These differences 
should be taken care of in the design and implementation 
of the data processing procedure to avoid errors. 

5. Data acquired from the 13.3 GHz system should be checked 
with respect to the phase shift error. Corrections to 
the variations of the aircraft parameters should be made 
as precise as possible in the data processing. 
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THE CALIBRATION DATA 


The laboratory calibration procedure for the 13.3 GHz scattero- 
meter was described in detail by Reid (ref. 15) . Therefore, only 
the calibration results will be given in the following. The cal- 
ibrate data presented here were obtained during August of 1976 
and were used in all the calculations in this document. Updating 
of these data are necessary whenever a' new calibration is made 

I 

for a major mission. 

A.l THE CALIBRATE CONSTANT 

As discussed in Section 3, the constant K which relates the trans- 
mit and receive signal levels to those actually recorded at the 
receiver output has to be predetermined in the laboratory. The 
input signal to the receiver for the calibration derives from a 
broad band noise diode. Following the procedure given by Reid 
(ref, 15), K was determined to be 125.6 dB. Although the cal- 
ibration procedure is almost identical for the 13.3 GHz and 1.6 
GHz scatterometers , there is a basic difference in the input 
sources used. While the broad band noise is used 'for the cali- 
bration of the 13.3 GHz system, the input source for calibrating 
the 1.6 GHz system comes from a single side-band modulator. Be- 
cause of this difference, the K values should be used with care 
in the data processing. 

A. 2 VARIABLE GAIN RESPONSE TEST 

The reason for this test was to verify that the changes in the 
output signal correspond to the changes in the variable gain 
amplifier settings. The test results were given in Table A-1. 

The gain difference A G was taken between any one of the outputs 
at X2, X5 and XIO gain settings and the output at XI gain setting. 
It is clear that the outputs of both channels correctly respond 
to the changes of gain setting. For example, when the (voltage) 
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TABLE A-l. - THE RECEIVER OUTPUT RESPONSES TO THE VARIABLE 

GAIN AMPLIFIER SETTINGS 


Gain 

Setting 


putput (dB) 


XI 

X2 

X5 


Channel 1 


-9.5 

-3.5 

+4.5 


AG CdB) 


+ 6.0 

+14.0 


Channel 2 


Output (dB) 


- 8.8 

- 2.8 

+5.2 


AG (dB) 


+ 6/0 

+14.0 


XIO 


+10.5 


+ 20.0 


+ 11 . 2 


+ 20.0 













Amplifier gain is increased by 10 times / the output power is 
observed to increase by 20 dB as expected. 

A. 3 THE HOLLOFF FUNCTION . 

The measurements of the sea (water) and land filter response as 
a function of audio-frequency (rolloff functions) were similar 
to that of determining the constant K. The results of these meas- 
urements were summarized in Table A-2 for both land and sea fil- 
ters. It is clear from this table that, for either land or sea 
filters, the frequency responses of the two receiver output chan- 
nels are very similar. The average frequency responses of the 
two receiver channels for both land and sea filters were plotted 
in Figure 3-3. Only the sea rolloff function was used in the cal- 
culations of the receiver output signal levels, because the land 

* 1 

rolloff filters were found not to give enough attenuation at low 
frequencies and were of little practical use. 


A. 4 THE RECEIVER NOISE FIGURE 

The measurements of the receiver noise figure were made at the 
frequencies of 2, 4, 6, 8 and 10 KHz. The results were summarized 
in the following table. 

TABLE A-3.- THE RECEIVER NOISE FIGURE 


Frequency 

(KHz) 

2 

4 

6 

8 

10 


Channel 1 
(dB) 

20.2 

17.7 ' 
16.2 

15.7 
15.7 


Channel 2 

21.7 

20.7 
18.2 
17.5 
17.5 


There is some differences in the noise figure between channel 1 
and channel 2 as noted from the table. 



TABLE A- - THE MEASURED RESPONSES (ROLLOFF FDNCTrONS) 
IN dB FOR THE 13.3 GHz SCATTEROMETER RECEIVER OUTPUTS 






gs 

tap 



Test 

Frequency 

(KHz) 

ChanneJLJL 

Channel_2 

Land 

Water 

Land 

Water 

10 

-0.2 

-0.3 

-0.2 

-0.5 

8 

-0.4 

-0.9 

-0.3 

-0.9 

4 

-0-6 

-3.1 

-0.5 

-3.5 

2 

-1.0 

-7.0 

-1.1 

-7.2 

1 

-1.6 

-12.5 

-1.3 

-12.6 

00 

• 

o 

-2.1 

-14.0 

-1.5 

-14.2 

0.4 

-4.1 

-20.0 

■-4.0 

-20.2 

CM 

• 

o 

-7.9 

-26.0 

-8.0 

-26.2 

0.1 

-12.9 

-32.0 

-13.0 

-32.2 

0.08 

-13.6 

-34.0 

-14.5 

-34.2 

0.04 • 

-19.6 

o 

• 

0 

1 

-20.5 

CM 

« 

O 

1 
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THE FERRITE MODULATOR 


TO derive the absolute value of the backscattering coefficient 
a® as a function of incident angle 6 from the backscattered 
Doppler signal return, it is necessary to maintain a reference 
calibrate signal, proportional to the transmitted power, in the 
receiver output. The 13.3 GHz scatterometer system uses a ferrite 
modulator to provide the reference signal required in the data 
reduction. The method is to amplitude modulate a sample of the 
transmitter power by a calibrate signal at 12 KHz and pass the 
sidebands of the modulation through the receiver network along 
with the backscattered signal. The ferrite modulator was installed 
in the channel which was not phase shifted. The single side-band 
power of the modulator as a function of the applied coil voltage 
was measured at several frequencies. The results were shown in 
Figure B-1. 

The ferrite modulator employs the Faraday rotation principle 
and operates as a variable attenuator attenuating the RF energy 
passing through the modulator (refs. 2 and 16) . If the maximum 
value of the rotation 6^ is sufficiently small, the output elec- 
tric field E^ from the 'modulator is given by (ref. 16). 

E = -ri [j^{6 ) + 2 J, (6_) sinw tl sinw^t (B-1) 

o2Lo m I m cj o 

where is the input electric field amplitude at frequency 

and Oi is the (calibrate signal) frequency at which -the modulator 
c ^ 

was driven by an audio oscillator. J and J, are the Bessel func- 

o 1 

tions. For small 6^^, << J^. Using different approach, Bradley 

(ref. 2) essentially arrived at the same relationship between the 
input and output voltages. 
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Figure B-1 -Ferrite modulator single sideband power 
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Thus the output voltage from the modulator can generally be 
v?ritten as: 


= a(l + ycosw^t) 

Substituting Eg. (3-16a) for V., Eq. (B-2) 

ak / 

^IF ~ 2 ~ jCf cos((w^ + 0)^) 

+ cos ((w^ ~ (U^) t + X^) + 

aki j 

= ^ Cf cos ({t«^ - ( 0 ^) t t X^) + 


(B-2) 

becomes 
t + X^) 

) 

— n — cosw^t , 

Jc jL o 

cos ( (w^ - w^) t + X^) 


+ ^ COSW t + c cosu t cosw t > (B-3) 

ki o kj^ o c j 

Two terms (product of calibrate and backscattered signals) inside 
the bracket in the above equation ^ 


YC_ cosu t cos ( (u t + w, ) t + x^) and yC cosu t cos((u ” u, ) t + X ) 

were omitted for the following reason. The ferrite modulator coil 
was generally operated at -.1,5 volts rms (S, C. Reid, personal 
communication) . From Figure B-1, the single sideband power was 
found to be -.*-90 dB below the power level of the local oscillator, 

X • 6 • f 


20 log ay ~ -90 dB (B-4) 

* 

Now the range of backscattered signal power at receiver input was 
from ~-80 dBm to -.-120 dBm depending on the incident angles. The 
product of the signal return power would be ~-170 dBm to ~-210 dBm 
which were way below the noise floor of --139 dBm given in Table 
2-3. Therefore, those two terms can be neglected and Eq. (B-3) 
represents the signal input to the mixer in channel 1. 
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THE MIXER MATHEMATICAL MODEL 
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THE MIXER MATHEMATICAL MODEL 


A mixer is essentially a nonlinear resistive element which com- 
bines the input signal with the local oscillator signal (the 
transmitter sample signal in the present case) and outputs sig- 
nals at frequencies harmonically related to the input and local 
oxcillator frequencies and their mixtures. The characteristics 
of the nonlinear resistance for a typical mixer is shown in 
Figure C-1. One of the terms in Eq. (B-3) can be regarded as 
the local oscillator signal; 

S^(t) = cos ^t (C-1) 

The remaining terms of Eq. (B-3) belong to S (t) . 

Now the power series representation of the current I^^ flowing 
in the mixer may be expressed in terms of the voltage across 
the mixer terminals (ref. 17) as; 

I- = a + a^ V + + a, . . . + a„ (C-2) 

Ira 0 12 3 n 


with 


V (t) = (t) + (t) 

The coefficient a generally decreases rapidly with n. Neglecting 
terms with n > 3, Eq. (C-2) becomes; . 

*■ , r 

By substituting the expressions for and S^, it is found that 
the signal given by Eq. (C-3) consists of DC terms as well as 

terms at the frequencies of ^o^ ^o ~^c' ^^o' 

2w, , 2w ±0),/ 2w ±co^r 2o)_ -2tD, / 2(o^ 2ta + u, 

d cd od oc o d o c o dc 

and 2 ( 0 ^ - scatterometer system under study, u)^ = 

13.3 GHz, {0 “12 KHz and Wj 5 9 KHz for aircraft speed < 200 knots. 

O Cl 
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During the subsequent amplifier stages all the terms at the fre- 
quencies > 12 KHz are filtered out. Therefore/ by retaining only 
the terms (except DC terms) at the frequencies < lo^/ the output 
signal voltage across the mixer load resistor R/ is: 

Vj^ = ■ ^ pCjP^V^oos + Xj) + o^c^pj^v^cos (a.^t - X^) 

+ cos ( (o)^ - 0 )^) t - X^) + cos ( (w^ + t + X^) 

+ a^C_YPiV cos ((to - to,) t - X ) + cos ((to + to.) t + X ) 

ci ± O Cu Si Cu Si 

2 ^a^YPi^V ^ \ 

+ a C-C^cos (2a,t + X^ - X ) t + — — cos t ( (C-4) 

r a a r a jc. c i 


The first/ second/ and last terms in Eq. (C-4) are the dominant 

terms. The power level of the local oscillator signal p^^V^costo^t 

is only about ~36 dB below the transmitted power (Figure 2-2 ) , 

which is much larger than the backscattered signal at ~460 m 

aircraft altitude. During calibration/ the power level of the 

calibrate signal at frequency to was set to be comparable to the 

c 

2 2 2 2 

backscattered signal/ i.e./U yPi V ~ a C-p,V or YPiV ~ C- 

1 o X 1 o 1 o X 

(or YP-|V ~ C^/ which also suggests that y<<D ♦ The ratio of the 

X O ci 

signal at frequency to that at to^ is: 


2“Cf0i''o " 2 


(C-5) 


and the ratio of ‘the signal at frequency 2to^ to that at to^ is: 


aCfP^Vo 


Pl^o 


Y«1 


(C-6) 
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Thus the signals at frequencies 2a)^ and can be neglected 

and Eq. (C-4) becomes: 

''im “ jcjCOS(a>^t + Xj) + C^oos - X^) 

2YPiV ) 

+ * ~ ^ ° cos(o t ( (C-7) 

ki c j 


with = a^a 2 R^. Following the same approach, it can be easily 
shown that the signal voltage output from the mixer in channel 2 
is: 


^2m = 


^2^2^2^o 


j-C^sin (w^t + X^) + C^sin (w^t - X^) j (C-8) 
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THE COMPUTER PROGRAM FLOW CHART AND LISTING 


In order to calculate the receiver signal output according to 
Eg. (3-26) , the incident angle 0^^, the distance R^^, and the area 
element a-s shown in Figure D-1, should be expressed in terms 

of the measured quantities such as along-track incident angle 
cross-track angle altitude h and the along-track cell length H. 

As mentioned before, the antenna measurements were made with a 

I ^ _ » 

conical cut around the X axis. In spherical coordinate system, X 
- axis becomes the polar axis. And the polar angle is 90° - d> and 
the aximuthal angle is 6^ measured in the YZ plane, in this coordi- 
nate system the center of the area element is located at 


X. = h secO”. tan4>. 
1 xp 1 

^i = h tan 6 . 

iP 

Z. = h 
X 


The Doppler ' s relation gives : 


sine 

X o 


2V 

X 


(Xf^ + + Z^^) ^ 


or 


sin0 = 


tan© . 


° (sec^0.„ tan^4>. + tan^0. + i) 




xp 


xp 


sinO. *= sin0 sec(|). 
ip ox 

The distance R^^, with the help of Eg. (D-1) , is given by 


(D-1) 


D-2 
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Figure D-1.— The geometry of an area element used in the numerical 

integration. 
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R. = h sec9. sec*. 
X xp ^x 


2 2 ¥ 
(cos’^<fi^ - sin 0^) 


(D-2) 


The angle of incidence 6^^ needed in obtaining a value for nvira- 
erical calculation is 


6. = cos ?r^ = cos (cos *. - sin 0^)^ 
X Rx ' ^x o' 


(D-3) 


The area element is given by 


AA. = AX. AY. 

X XX 


(D-4) 


with 


AX. 

X 


X.., X., “ h 

X2 xl 


sin(j> 


i2 


~t 2T . 2„ 

(cos (^^2 ~ ®o^ 


2 

(cos - sin0^) 


(D-5) 


£cos0 


and 


AY. = 


4 v.2„ ^ 


^ (cos'^d). “ sin'^O^) 

' ^x o 


(D-6) 


The angles and 4>j^2 related to 4)^^ (all units in degrees) by 


4*^0^ — 4>^ 0.5 


(D-7) 


4>i2 = 4>i + 0.5 ' (D~8) 

The above expression for R. , 0., and AA. were used -in the computer 
program for calculating the receiver signal output. 


The flow charts of the main program and the subroutine GAAT are 
displayed in Figures D-2 and D-3 respectively. The listing of 
the whole program is given in Table D-1. Three subroutines RINIT, 
RWRITE/ and RREAD are used to initialize, to write data into, and 


D-4. 
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START 


ENTER CONSTANT PARAMETERS 


COMPUTE AREA, 
DISTANCE, AND SIGNAL 
CONTRIBUTION FROM THE 
AREA 


READ AND STORE ANTENNA 
PATTERN DATA 


ENTER THE DO LOOP OF 
SIGNAL COMPUTATION 
AT SELECTED INCIDENT 
ANGLES 


CONTRIBUTION \ 
FROM ALL CROSS-TRACK 
AREA ELEMENTS 
ACCOUNTED FOR? ^ 


CALCULATE DOPPLER 
FREQUENCY AND BANDWIDTH 


COMPUTE AND STORE TOTAL 
OUTPUT SIGNAL FOR AN 
INCIDENT ANGLE 


ENTER THE DO LOOP OF 
COMPUTING CONTRIBUTION 
FROM VARIOUS AREA ELEMENTS 
IN THE CROSS-TRACK 
DIRECTION 


^ COMPUTATION 
FOR ALL SELECTED 
INCIDENT ANGLES 
. COMPLETED? 


CALL GAAT 


DERIVE ANTENNA GAIN 



OUTPUT RESULTS 



Figure D-2.— The MAIN program flow chart. 










Figure D-3.— Subroutine GAAT flow chart. 
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TABLE D-l.-THE COMPUTER PROGRAM LISTING 


1« 

2 *. - 

3« 

5« 

7*- - 
8# 


9« 

1 C>!‘ 
11 *. 
12 * 
13* 


14* 

15 *- 

16* 

17* 

18* 


19* 

2f* 

21 * 


2 2 * 
-2 3*- 


24* 

25* 

26* 

-2^* 


28* 
29* 
3 0* 
3-K 


3 2* 
3 3* 
34* 
35fr. 


36* 

37* 

38* 

39* 


4 0* 

4 1# 

4 2* 

43* 

44* 

45* 

46* 

-4-7*- 

48* 

4 9* 

5 0* 
51* 
5 2* 
5 3* 
54* 
- 55 * 


56* 
57* 
58* 
59* - 
6 0 * 
61* 
62* 


niM£ NSION theta ( 8 ) , P CW ( 1 6) ,R0FF (8 ) »BUF ( 181 ) 

qaTA the TA/5r , 10, , 1 5 .|2n, , 30 . ?4 L. ,5u . » 6C»/ 

DATA R0FF/I6.fa,ir. ,8,uR,6,t},4,D,2,9,2,2»l,8/ 
rOMHON/B/IPACD, IX ,I Y ,1XMAX, IY*A X 
C ENTER CONSTANT PAPAPETEPS. 

PIE;3.1416 

VAxtol” 

SLAND=0. 02255 
rELEN=25. 

PI=PIE/1&0, 

SCAL=125.6 
C GET DISC AREA 

TAIL RINIT ( IDADD, NWR ES ) 

K-=IDADO 

PEAL ( 1 )I YMAX,IXMAX, D El , DB2 , DPMAX 
IX=( IXMAX*?. J-H 
TYr( lYMAX*? , ) +l 

C— PUT-ANTENNA GAINS -OM-OISC-* 

no 100 I=1,IY 

RE AD < i )N , Blip 

WRITE ( 6, 4u^0) (BUF (K J ,Kr 1 ,IX ) 

CALL RWPITE(K,BUF ,IX ,NSTAT~) 

5 IF (NSTAT ,E0. 1 )60 TO 5 

IF (NSTAT.NE. JJWRITE ( E,TG00)N 

t/ — i/x T Y 

— lO n 00 N T I N UE- 

C ENTER MAIN LOOP OF COMPUTATION. 

C POWER COMPUTATION FOR EACH INCIDENT ANGLE, 
no 5''U 1 = 1,8 

TH TA =TFtE-T-A-<-I-> *PI 

F0 0P=2.*VA=^SIK(THTA )/SLAMD 

DF=2 .*VA*{ COS(THTA) *>!3 ) *CEL EN/H A/ SLA HD 

CLAMO=SL AMD**2 

CPI= ( 4 ,*PIF-1**3 


SUPFrsJ. 
SUMA =0. 


C CROSS TRACK COMPUTATION. 

DO 4 CQ-J = lr 7 

PHI = ( J-n*PI 
PHIL=( ( J- 1 ) - 0 . 5 )*PI 
PHIR=( ( J- 1 ) + 0 , 5 )*PI 

F-AC=SQRT ( (COS (PHI P) * »2 ) ~(S-I N ( 7 H T A-) ** 2 -H 

XPHR=SIN(PHIRJ/FAC 
XPHL=SIN ( PHIL) /FAC 
PX=HA*( XPH 9 -XPHL ) 

nv=CELFN*COS( THTA >/-F AC 

DA=LX*DY 

R=HA/FAC -r 

RRrRi:t« 4 , / ■ 

TH 1 NC=AC 0 SIFAC) 

THINC=THINC/Pi 

SIGHA = { 163 . 2 /(THINC*>!C. 76 ) )- 30 . 

AI lA SIN ( SIN { ThT A ) /CO S(PHI ) ) 

A 4 =A I /PI 

CALL GAAT(FGL,FGR,AGL,AGR,J,AI) 
SIG= 10 .**I SIGMA/n. ) 

SUMA=SUMA + ( { AGL-^ AGR ) <S I 6*0 A /PR » 

SUMf rsUMF-M (FbL + FGR ) >!S T G*D A /PR >— 


WRITE ( 6, 6CnO)FGL , FGR ,AGL , A 6R , SU HF , SU MA , SIGMA, OA , RR 
40r CONTINUE 

C CALCULATL POWER OUTPUT FOR AN INCIDENT ANGLE. 
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TABLE D-l.“ Continued 


63« 
64# 
65# 
66 # 
67# 

•6 8 #" 
69# 

7 0# 
71# 
-7?#- 

7 3# 
74# 
75# 

.76#* 
77# 
78# 
79# 

-8 0 # 
81# 

8 2 # 

8 3# 
84# - 


FPAC =SUMF#CL AMD/CPI/CF 
- APAC =SUMA#CLAHD/CPI/CF 
FPAC=in«#AL0G10(FPAC ) 

APAC=lC.*ALOGl-?( APAC 1 

WRI7E(6,6 7nn)FPAC,AP AC,FD0P,DF, CLAHD fCP I,SCAL,DX,DY 
— POW( I >:FPAC“OOFF-< I) +-SCAL- 


500 


P0W(I + 8)=APAC-R0FF{I J-t-SCAL 
C0^aINUE 
WRITE( 6» inOD) 

WRITE! 6 , 20nc MPOW M- >-<!-= 1 1 8 -J 

WRITE! 6f 3rjno) 

WRITE! 6, 2t3''-'D) !P0W1I> ,1 = 9,16) 

STOP 

■lOOU— -FORMAT! 2M- r’-INCI-0£N T- ANGLE-AT ’ ,-3X"r-^6-J-rHX-rr*-l-0'VlCX-,"-*-l-S* r 
ll'^X,*2u», ICX,’ 3D*, 10 X ,*4 0 MUX, ’60 Si OX, *6u'S//) 

2000 FORMAT! 16X, 8012,41 

3UC^ F0RHAT(2H , ’INCIDENT ANGLE AT * , 2X , *- 5 ’ , !DX , 10 * , 9X , 15 * , 

1.0 V- • - xn 9X , J-4 0 9 X , 50 S-9 X , 


-1-9X-, *-20* , 94(r^*--50» 
40CC F0RMAT(20F6, 1) 
600J F0RNATt9D12.4) 
70CU F0RMAT!2H ,’TAPE 
• E-ND- -- 


PROBLEM AT ANGLE*, 15) 


1 # 

2 # 

3# 

— 4#- 


5# 

6 # 

7# 

9# 

10# 

U# 

• 12#" 
13# 
14# 
15# 
16# 
17#‘ 

le# 

19# 
2 0 # 
21#- 
22 # 

24# 

• 2S#- 
26# 
27# 
28# 
2 ?#- 
3 0# 


SUBROUTINE 6AAT!FGL » FG» , AGL ,AGR , J»AI ) 

g?SrWE§!!Pl{?jn;iS3fE*!5^!}^6uri(,,n,rBur2,ui. 

M.J NE - T H E-P 0 SH 1 0 N S' ■ O^F">T-F! E--A N T-E-N N A-f 0 0 T-P RI N-T"S 

IRSIXMAX 
ri =T XMAX4R 


an “a /> • I « 

1L=IXMAX42 
TM=IXMAX+1 
-I A =AlWr-!‘AI-4- 
1APPSIA+ 1 

tel *«* teV ,111 1 I 




..lUlllfJA- ■* 


TAFl=iyMAX+IA 
TAA1=I YMAX-I A 

A F 2 s I AF 1# 1 — 

IAA2SIAA1-1 
lAPlSl AA 1#T X + IDADD 
IAP2 = 1 AA2#T X*IDAOO 
IFP1SIAF1#IX+I0A0D 
TFP2 = IAF2#IX + IOAOD 

-• ir (1 ST InE I 0 1 WRI TE ! 6 ♦ pOQ > 1 A— ' 

CALL RR£AD(IAP2, ABUF 2,IX,1ST) 

lF(ilTlN?tMwRIT?!6, 1000) I API 
••— ^A U RRE AD ! I FP 1 ,F BUF‘ 1 , 1 X rIST-1 
TFIIST.EQ.DGO TO 7 
TF nST.NE.P)WRXT£ (6, 1U^0)1A 
CALL RREA0!IFP2,F8UF2,XX,IST ) 
■“•••IF ! 2 ST ,£0 , 1 )60 ■ TO 8* - 

XF (1ST. NE.nj WRITE ! 6, 1000)1 API 
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TABLE D-1 


Continued 


31« 

3 2* 

■ 33* -- 
34* 

35* 

36* 

. 37* 

38* 

39* 

4 0* 

- 41* 

i» 2* 

43* 

44* 

45<c 

46* 

47* 

48* 

-49* 

5G* 

51* 

5 2* 

S3* 

54* 

5 5* 

56* 

-5-7* 

58* 

59* 

60* 

- 61 * 

6 2 * 

6 3* 

64* 

-- 65 * 

66* 

67v 

68 * 

-69^ 


WRITE (6, 20"0nABUFl ( ,tCl 
WRITE(6,?C'?rHABUF2(K),K=i 
— -WRITE<6,2C1''f'MFBUFl ( 4),K=1 
WRITE ( 6, 2d'’C-) (FBUF2 ( H) , Kri 
OBTAIN THE ANTENNA GAIN VALUES 
IFtJ.EO.lJGO TO 10 

BEAM antenna GAIN 

APL1=ABUF1 ( IL-J) 

APL23ABUF2 t IL-J ) 
APR1:ABUF1(IR+J) 
APR2=ABUF-2{I-P+J->- 


IX 1 
IX ) 
IV >-■ 
IX ) 


_C -A F-T 


C FORE 


APL1L=10.**( APLi/in 
APL2L=ia.«*(APL2/lC 
APPiL=10.**{ APRl/ in 
-APRi,L = l&.**« AFR2/ 1C ■ 
BEAM ANTENNA GAIN 
FPLl^FGUF 1 (IL-J ) 
FPL2=FEUF2(IL-J) 
-FPR1=FFUP-1 ( IR+J J 


) 

J 

) 

- 1 - 


FPRirFBUF 2 ( IR+J ) 

FPLIL=1G.**(FPL 17 10. 1 
FPL2L= IG.** ( FPL 2/ 12 . ) 

FPRlL = in.v*(FPRi/ I-'- .-> 

FPR2L=10.*’"(FPR2/10. ) 

LINEAR INTERPOLATION. 

AGL={((AI-IA>/(IAPP-IA) ) *( A PL2L -APLl L ) ) +APL IL 

A6P= ( ( ( AI-I A ) / ( I APP— lA > )*(-A PP?L— APRIL-} )--»AP-R-lL— 

FGL={((AI-IA)/(IAPP-]A1 ) *( FPL2L -FPL 1 L ) ) +FPL IL 
F6R={((AI-IA)/(IAPP-IA})*(FPR2L-FPR1L))+FPR1L 
GO TO 15 
CONTINUE- 


-Id- 


C AFT BEAM ANTENNA GAIN, 
AP1=ABUF KIM) 
AP2=ABUF2( IM> 

A p 1 r 1 0 , * =fr.GA P-l-/-l-iI,-> 

AP2=10.**( AP2/1C. ) 
C FORE BEAM ANTENNA GAIN, 
FPI-FBUF KIM) 
-FP2=FeUF?( IH ) 


ALONG TRACK 


ALONG TRACK 


7C* 

71* 

72* 

73* 

74* 

75* 

76* 

77* 

78* 

C LINf 
15 

79# 


8 0* 

10 no 

81* 

20ru 

82* 

4QQC 

8 3* 



FPUl J.««(FP1/10. ) 

FP2=10.**(FP2/13. ) 

iR INTERPOLATION. 

A6L= K ( AI -I A )/ (1 APP- lA) J * ( AP2-AP1 ) ) + API 
AGR=D. 


F6L=(< (AI-IA)/(IAPP-1A» )*{FP2-FP1)) + FP1 
F 6 R- 0 . 

...WRITE { 6, 40CO) AGL♦-AGR.,F6L-^f.6R.^AP-U-AP^-,FJ>.♦itKp-2^-J^-I-r-IA^•t•d■ 
PE TURN ' 

rORHATtlH ,*DISC RE A C PROBLEM FOR ANGLE*, 15) 
F0RMAT(2QF6.1) 

FORMATtlH ,9F&.2,2I1C> 

END . 
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to read data out of the high speed drum storage. These subroutines 
are stored in the local computer program library and explained 
in the JSC CAD Procedure Manual, Exec 8 System General User Infor- 
mation . 


origwaIj,Mgek 
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